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STUDIES  OF  PLASMA  INSTABILITY  PROCESSES  EXCITED  BY  GROUND 
BASED  HIGH  POWER  HF  ('  HEATING")  FACILITIES 


The  work  is  driven  by  the  intersection  of  two  lines  of  activity,  one  experimental,  the 
other  theoretical,  that  converge  on  an  importantly  new  view  of  the  physics  of  how  high 
power  HF  radio  waves  interact  with  collisional  plasmas,  such  as  the  earth's  ionosphere. 

Its  goal  is  to  seize  this  opportunity  to  amalgamate  the  complementary  US  and  Russian 
recent  findings,  to  strengthen  the  foundation  of  plasma  physics  underlying  future  work  in 
the  growing  field  of  ionospheric  HF  heating. 

The  main  new  results  can  be  summarized  as  following: 

1.  Two  sets  of  observations  of  suprathermal  electrons,  produced  by  the  interaction  of 
powerful  radio  wave  with  the  ionosphere,  as  seen  using  the  incoherent  scatter  radar 
(ISR)  technique  are  presented.  The  observational  data  are  compared  with  the  theory  of 
multiple  acceleration  of  electrons  in  the  strongly  excited  resonance  region  near  the 
reflection  point  of  the  powerful  radio  wave.  The  structure  of  the  wide  perturbed  region 
filled  with  energetic  (10-20)  eV  electrons  is  determined.  The  size  of  this  region  along 
the  Earth  magnetic  field  is  shown  to  be  100  km.  The  full  power  going  to  the 
accelerated  particles  is  determined:  it  is  6-8  kW  or  (4-6)%  of  the  entire  HF  radiated 
power.  The  power  carried  by  the  suprathermal  electrons  flux  escaping  into 
magnetosphere  is  of  the  order  1  kW. 

2.  The  of  electrons  near  the  O-wave  reflection  point  by  combimed  action  of  upper  hybrid 
and  Langmuir  turbulence  is  considered.  It  is  demonstrated  that  the  calculations  of 
electron  acceleration  in  Langmiur  resonance  layer  without  electron  collisions  lead  to 
controversial  result.  Kinetic  theory  of  multiple  acceleration  of  suprathermal  electrons 
is  developed.  The  energy  losses  of  electrons  on  the  excitation  of  molecular  vibration 
level  are  taken  into  account.  The  model  calculations  describing  the  general 
peculiarities  of  the  suprathermal  electron  space  energy  distribution  are  performed. 

3.  Irregularities  of  electron  temperature  and  plasma  density  in  the  F-region  of  the 
ionosphere  are  strongly  elongated  magnetic  field  due  to  dramatic  difference  in 
transport  coefficients  along  the  magnetic  field  and  transverse  to  it.  Their  scales  along 
the  magnetic  field  determined  previously  neglecting  plasma  drift  effects  are  rather  long 
reaching  ten  or  even  several  tens  kilometers.  Such  elongated  irregularities  should  be 
affected  by  ionospheric  drifts  perpendicular  to  magnetic  field,  but  up  to  date  there  was 
no  theoretical  estimations  for  this  influence.  Two  main  problems  are  solved:  constant 
drift  effects  on  irregularities  caused  by  the  sources  fixed  in  the  ionosphere  (for 
example,  narrow  radio  beam),  and  gradient  drift  effects  on  irregularities  created  by  the 
source  moving  with  plasma. 
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1.  Experimental  Studies  of  Energetic  Electrons  in  Ionosphere  Using  ISR 
Technique  and  Comparison  of  Observational  Data  with  Multiple  Acceleration 

Theory 


Abstract 

Two  sets  of  observations  of  suprathermal  electrons,  produced  by  the  interaction  of 
powerful  radio  wave  with  the  ionosphere,  as  seen  using  the  incoherent  scatter  radar  (ISR) 
technique  are  presented.  The  observational  data  are  compared  with  the  theory  of  multiple 
acceleration  of  electrons  in  the  strongly  excited  resonance  region  near  the  reflection  point 
of  the  powerful  radio  wave.  The  structure  of  the  wide  perturbed  region  filled  with 
energetic  (10-20)  eV  electrons  is  determined.  The  size  of  this  region  along  the  Earth 
magnetic  field  is  shown  to  be  100  km.  The  full  power  going  to  the  accelerated  particles  is 
determined:  it  is  6-8  kW  or  (4-6)%  of  the  entire  HF  radiated  power.  The  power  carried  by 
the  suprathermal  electrons  flux  escaping  into  magnetosphere  is  of  the  order  1  kW. 

1.1  Introduction 

Observations  of  fast  suprathermal  electrons,  due  to  the  interaction  of  powerful 
radio  waves  with  the  ionosphere,  using  incoherent  scatter  radar  (ISR)  technique  were  first 
proposed  and  realized  by  Carlson,  Wickwar  and  Mantas  in  1972  (see  Carlson  et  al. 
[1982]).  The  powerful  O-mode  wave  in  the  resonance  region  near  the  reflection  point 
effectively  excites  plasma  waves  -  natural  oscillations  of  the  ionosphere  plasma.  Plasma 
waves  due  to  the  nonlinear  cavitation  process  accelerate  suprathermal  electrons  to  energies 
in  excess  of  e  ~  1  0+20  eV,  i.e.  two  orders  of  magnitude  higher  than  the  thermal  electron 
energy  Te  *  0.1  eV.  The  fast  electrons  propagate  in  the  ionosphere  to  large  distances,  10- 
100  km  from  the  acceleration  region.  Those  electrons  excite  nonthermal  plasma  waves, 
which  were  detected  by  ISR  [Carlson  et  al.,  1982]. 

We  discuss  here  the  results  of  two  sets  of  experiments.  The  experimental  data  are 
here  compared  in  detail  with  theory  based  on  the  mechanism  of  multiple  electron 
acceleration. 

Carlson  et  al.  [1982]  gave  a  qualitative  explanation  of  their  observations.  They 
pointed  out  the  significant  role  played  by  multiple  scattering  of  electrons  in  the  neutral 
atmosphere.  The  quantitative  theory  of  multiple  acceleration  was  proposed  by  Gurevich  et 
al.  [1985].  The  comparison  of  the  theory  and  observations  can  help  reach  a  still  higher 
level  in  our  understanding  of  a  number  of  space  plasma  phenomena  in  which  accelerated 
electrons  play  an  important  role. 


1 .2  Observations 

The  most  sensitive  and  direct  ground  based  means  of  detecting  10+20  eV 
suprathermal  electrons  is  by  incoherent  scatter  plasma  line  observations  of  the  Langmuir 
waves  produced  in  the  ionosphere  by  these  electrons  [Perkins  and  Salpeter,  1965].  This 
technique  has  been  used  to  study  photoelectron  fluxes  in  the  ionosphere  [Yngvesson  and 
Perkins,  1968].  Carlson  et  al.  [1982]  applied  the  ISR  technique  to  measurements  of 
suprathermal  electrons  produced  by  the  interaction  of  powerful  radio  wave  with  the 
ionosphere  at  Arecibo.  We  present  here  two  types  of  experimental  data.  The  primary  and 
secondary  data  were  gathered  on  the  nights  20  May  1972  and  13  July  1992  respectively. 
The  primary  data  were  published  by  Carlson  et  al.  [1982],  here  we  will  revisit  those  results 
and  supplement  them  with  the  complimentary  set  of  secondary  1992  data  [Djuth  et  al., 
1996], 
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For  the  1972  Arecibo  data  here,  a  HF  dipole  feed  hung  over  the  1000-foot  diameter 
ISR  radar  dishes,  for  which  the  138  kW  transmitter  at  7.63  MHz  delivered  about  30  W/m2 
to  the  HF  reflection  height.  For  the  1992  Arecibo  data  here,  a  separate  48  array  of  log 
periodic  antennas  17  km  away  from  the  ISR  radar  dish,  for  which  the  300  kW  transmitters 
at  5.1  MHz,  into  50%  efficiency  transmission  lines,  delivered  about  40  W/m2  to  the  HF 
reflection  height. 

For  20  May  1972,  strong  plasma  waves  were  excited  near  the  reflection  point  at  the 
altitude  of  z  =  285  km.  Night  time  plasma  line  intensities  were  observed  to  be  enhanced  by 
a  factor  of  10-100,  over  a  range  of  altitudes  which  extended  to  below  250  km.  They 
relaxed  to  their  normal  level  within  the  electron  transport  time  after  the  transmitter  was 
turned  off. 

The  430  MHz  incoherent  scatter  radar  was  used  to  diagnose  the  ionosphere,  which 
allowed  measurement  of  altitude  profiles  of  the  background  plasma  temperature  Te  and 
electron  concentration  N.  The  plasma  line  echo  from  the  plasma  waves  excited  by  the 
suprathermal  electrons  provided  the  altitude  of  a  set  of  plasma  frequencies  fp  between  5 
and  7.5  MHz,  with  an  altitude  resolution  of  1.5  km. 

For  a  given  diagnostic  radar  wavelength  Ar  the  main  scattered  signal  comes  from 
ionosphere  plasma  waves,  whose  wave  vector  is  directed  toward  the  radar,  and  the  phase 
velocity  is  equal  to  the  velocity  of  suprathermal  electrons  v 

v„=±a,4  a) 

Here  fp  =  (e2N/47tm)12  is  the  local  plasma  frequency.  For  the  Arecibo  radar  Xr  =  70 
cm  and  from  eq.  (1)  one  can  obtain  that  the  energy  of  electrons  is 

*  =  O.354V  (2) 

The  intensities  of  plasma  line  echoes  were  shown  in  Fig.l  of  Carlson  et  al.  [1982], 
A  significant  enhancement  of  radar  echo  intensities  is  evident  at  large  distance  below  the 
reflection  point  of  the  radio  wave  at  285  km.  The  upgoing  one-dimensional  fluxes  of 
suprathermal  electrons  at  the  heights  266  km  and  256  km  determined  from  these  data  are 
shown  in  Fig.3  of  Carlson  et  al.  [1982]. 

The  additional  data  were  gathered  on  the  night  13  July  1992.  The  reflection  altitude 
was  295  km,  with  strong  enhancements  of  plasma  line  intensity  observed  in  the  vicinity  of 
the  reflection  level.  An  example  of  plasma  line  echoes  collected  at  the  heights  337-367 
km,  where  the  ISR  beam  intersected  the  flux  of  the  accelerated  electrons,  is  shown  in 
Fig.  1.1.  The  echoes  were  obtained  both  below  and  above  the  maximum  of  F-layer. 
Electron  density  Nmax~  5.5- 105  cm'3  was  reached  at  F-maximum  heights  z  =  350  km.  The 
observed  plasma  line  echoes  correspond  to  upgoing  electron  flux.  Note  particularly  here  a 
significant  enhancement  of  plasma  line  echoes  even  at  altitudes  more  than  70  km  above 
the  reflection  layer. 


1 .3  Brief  Outline  of  the  Theory 

In  ionosphere  modification  experiments,  electrons  gain  energy  in  a  strongly 
disturbed  Langmuir  resonance  layer  near  the  reflection  point  of  the  O-mode  wave.  The 
essential  point  of  multiple  acceleration  theory  is  a  large  fraction  of  fast  electrons,  after 
leaving  the  acceleration  layer  will  return  back  due  to  collisions  with  neutral  particles,  and 
thus  can  gain  additional  energy.  The  process  can  be  repeated  many  times.  Thereby,  a  wide 
region  around  the  acceleration  layer,  which  is  elongated  along  the  geomagnetic  field,  can 
be  filled  with  strongly  heated  magnetized  suprathermal  electrons.  The  elongation  along  the 
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Earth's  magnetic  field  can  reach  an  order  of  a  hundred  kilometers  because  it  is  determined 
by  two  factors:  the  large  mean  free  path  of  the  suprathermal  electrons  and  a  small  part  of 
electron  energy  is  lost  in  one  collision.  The  extent  of  upward  elongation  will  be 
considerably  greater  than  in  the  downward  direction,  because  the  mean  free  path  above  the 
reflection  height  is  considerably  greater  than  below  it.  Consequently,  a  significant  flux  of 
suprathermal  electrons  moves  to  the  magnetosphere,  and  can  even  reach  the  magnetic 
conjugate  point. 

The  acceleration  process  due  to  the  multiple  crossing  of  acceleration  layer  is 
averaged  and  in  final  form,  depends  on  two  scalar  factors  only:  the  full  power  density  P 
absorbed  by  fast  electrons  in  the  acceleration  layer,  and  a  characteristic  parameter 
describing  the  effectiveness  of  the  acceleration  inside  the  layer  Tef.  They  are  related  to  the 
details  of  the  accelerating  process,  effective  number  of  cavitons,  their  width,  and  so  on. 

The  distribution  fiinction  of  suprathermal  electrons  in  the  theory  of  multiple 
electron  acceleration  can  be  presented  in  the  simple  form  [Gurevich  et  al.,  1985] 

dz  "  (3) 


f0(E,z)  =  CK0 


exp 


'Zrcosa 


Where  Ko  is  a  modified  Bessel  fiinction,  s  is  the  electron  energy  and  Tef  is  the 
effective  temperature  of  suprathermal  electrons,  while  a  is  the  angle  between  the  vertical 
and  geomagnetic  field.  The  normalization  constant  C  is  directly  proportional  to  the  power 
density  P  of  the  HF  wave  absorbed  by  the  suprathermal  electrons: 

C- - tP  (4) 


4  nXf 


(S/Y2 


Here  5  is  the  average  fraction  of  electron  energy  lost  in  a  single  collision  with 
neutral  molecules. 

As  given  in  eq.  (3)  the  acceleration  layer  is  assumed  to  be  located  at  z  =  0,  while 
the  Lc±  factor  is  the  characteristic  relaxation  length  of  suprathermal  electrons  in  the  upward 
(+)  and  downward  directions  (-) 

Z(z)=\k<’„¥PsY,\  (5) 

Where  Nm*  =  Nm±(z)  is  the  neutral  density  above  (+)  and  below  (-)  the  layer,  atr(e) 
is  the  total  transport  cross-section  of  electron-neutral  collisions,  and  8=  am/  <jtr,  where  o/„ 
is  the  total  cross-section  of  inelastic  collisions,  which  includes  ionization  by  electron 
impact.  In  <j,r  and  ain  the  collisions  with  all  neutral  components  are  considered. 


1 .4  Comparison  of  the  Observations  with  Theory 

In  the  experiment  on  20  May  1972,  upgoing  (in  the  radar  look  direction)  fluxes  of 
suprathermal  electrons  in  the  energy  range  10  -5-  17  eV  at  the  heights  z/  =  256  km  and  Z2  = 
266  km  were  determined.  The  reflection  point  was  at  zo  =  285  km.  Comparison  with  the 
theory  requires  determination  of  the  electron  flux  along  radar  direction  er: 

=  (v,er)/(v) 

We  introduce  net  upgoing  and  downgoing  fluxes  by  integrating  the  flux  Jr  over 
corresponding  angles  in  velocity  space.  Taking  the  angle  between  the  geomagnetic  field 
and  vertical  radar  ray  in  Arecibo  as  a  =  40°  we  obtain  the  upgoing  J+  and  downgoing  J 
fluxes  in  energy  interval  de. 
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The  results  of  calculations  are  presented  in  Fig.  1.2  for  the  given  heights  and 
different  Tef.  It  is  apparent  from  the  figure,  that  in  the  energy  range  10  h-  20  eV  reveal  a 
rather  flat  spectrum.  Besides,  there  is  no  strong  difference  between  the  fluxes  J  at  different 
temperatures  Tef. 

On  the  other  hand,  as  follows  from  calculations  the  fluxes  are  effectively 
diminishing  with  distance  from  the  acceleration  layer.  This  is  because  electron  energy  is 
lost  in  inelastic  collisions  with  the  neutral.  Equation  (7)  implies  that  the  fluxes  of  upgoing 
electrons  are  smaller  than  the  downgoing  fluxes  below  the  acceleration  region.  Their 
relation  depends  upon  the  angle  a  between  radar  ray  direction  and  the  geomagnetic  field  - 
for  Arecibo  at  a  =  40°  J+/J.  w  1.5.  Comparison  between  the  theory  and  observations 
presented  in  Fig.  1.3  shows  a  reasonable  agreement  between  those  two.  In  fact,  the 
behavior  of  the  spectrum  of  the  electron  flux  at  different  heights  is  consistent  with  the 
theory  for  any  of  2V  applied.  Taking  into  account  the  absolute  values  of  upgoing  flux  J+  = 
(4-f-8)xl0  el/cm  s  eV  at  characteristic  energies  1  0h-1  5  eV  obtained  by  Carlson  et  al. 
[1982],  one  can  find  the  absorbed  power  Ws  of  the  HF  wave  converted  into  acceleration  of 
the  suprathermal  electrons 


for 

Tef=  5  eV 

Ws~5.6kW 

for 

Tef—  7.5  eV 

Ws~6.2  kW 

for 

Tef=  10  eV 

Ws~8.5  kW 

In  these  calculations  the  Arecibo  heater  beam  was  taken  roughly  103  km2  for  a  full 
absorbed  power  Ws  =  PS.  The  factor  P  was  determined  from  calculations  of  the  electron 
flux  made  at  different  distances  from  the  source,  and  from  the  results  presented  in  Fig.  1.3. 
It  is  apparent  that  the  dependence  of  Ws  on  Tef  is  not  very  strong.  Note  that  eq.  (7)  allows 
us  to  estimate  the  upgoing  electron  flux  at  z  >  500  km  which  escapes  to  the  magnetosphere 
due  to  the  absence  of  scattering  collisions. 

Figure  1.4  shows  the  height  distribution  of  plasma  line  intensity  measured  during 
the  second  experiment,  which  is  compared  with  the  height  dependence  of  the  distribution 
function  of  suprathermal  electrons  fo(e(z),z).  Here  energy  e(z)  is  determined  through  the 
measured  electron  density  distribution  (Fig.  1.1)  by  using  eq.  (2).  One  can  see  a  sufficient 
agreement  between  the  theory  and  observations.  Note,  that  the  normalized  height 
dependence  of  the  distribution  function  practically  does  not  depend  on  the  parameter  Tef 
(see  Fig.  1.4). 

1 .5  Discussion  and  Conclusions 

It  was  shown  that  the  ISR  plasma  line  measurements  are  in  reasonable  agreement 
with  the  theory  of  multiple  electron  acceleration.  From  the  observational  data  and  their 
comparison  with  the  theory  it  follows  that:  In  experiments  where  powerful  radio  waves 
interact  with  the  ionosphere  a  large  number  of  suprathermal  electrons  are  generated  in  the 
energy  range  up  to  20  eV. 

The  suprathermal  electrons  are  observed  over  a  wide  altitude  region  of  the  order  of 
a  few  tens  km  both  below  and  above  the  acceleration  layer,  centered  near  the  reflection 
altitude  of  the  powerful  radio  wave.  While  not  observed,  they  must  extend  well  above  this 
altitude,  into  the  plasmasphere/magnetosphere,  and  even  to  some  degree  into  the 
magnetically  conjugate  region. 

In  power,  Ws  ~  5.6  h-  8.5  kW  (4  h-  6)%  of  the  full  transmitted  heater  power  goes 
into  the  acceleration  of  suprathermal  electrons.  This  energy  is  dissipated  due  to  the 
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generation  of  optical  emissions,  ionization  and  heating  of  ionized  and  neutral  components 
of  the  ionosphere  plasma.  The  characteristic  dissipation  length  for  suprathermal  electrons 
depends  strongly  on  the  altitude.  For  z  =  250  km  it  is  about  10  km,  while  for  z  ~  300  km  it 
is  about  30  km. 

The  flux  of  suprathermal  electrons  into  the  plasmasphere/magnetosphere  depends 
strongly  on  the  height  zq  of  the  acceleration  layer.  For  the  studied  cases  zq  =  285  km,  and  zo 
=  295  km,  and  Ws  =8.5  kW,  Te/  =  10  eV  the  flux  into  magnetosphere  in  the  energy  range 
10  4-  20  eV  is  about  0.9  kW.  This  estimate  is  consistent  with  earlier  calculations  of 
photoelectron  energy  and  flux  escape  into  and  loss  to  the  plasmasphere/magnetosphere 
[Mantas  et  al.,  1978], 

We  conclude  that  plasma  line  observations  combined  with  the  theory  enables  us  to 
obtain  significant  information  about  the  acceleration  of  the  suprathermal  electrons  in 
ionosphere  modification  experiments.  For  the  first  time  the  structure  and  size  of  the 
perturbed  region  filled  with  suprathermal  electrons  along  with  the  full  power  going  into  the 
accelerated  electrons  is  determined  by  a  theoretical  calculation  consistent  with  an 
observational  data  set. 

Note  that  we  had  the  possibility  to  use  here  only  a  small  part  of  ISR  data  on 
ionosphere  existing  from  experiment  1992.  The  further  elaboration  of  this  data  and  their 
comparison  with  the  theory  is  of  a  significant  interest  and  will  be  definitely  quite  fruitful. 

We  suggest  also  these  experiments  be  repeated  to  obtain  more  information  about 
the  main  features  of  acceleration,  and  its  dependence  upon  the  reflection  height  and 
ionosphere  conditions.  Note  that  such  observations  have  as  yet  been  done  only  at  Arecibo, 
though  we  expect  that  the  latitude  dependence  of  the  effect  could  be  very  significant. 

At  the  same  time,  the  agreement  that  found  between  the  theory  and  observations 
make  it  possible  to  plan  more  detailed  future  experiments.  These  will  thus  help  to  reach  a 
much  better  understanding  of  the  physical  mechanisms  of  electron  acceleration  and 
nonlinear  processes  in  ionosphere  plasmas. 
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Figure  Captions 

Figure  1.1  Spectra  of  the  downshifted  plasma  line  measured  at  Arecibo  13  July 
1992  at  8:30  LT. 

Figure  1.2.  Upgoing  flux  of  suprathermal  electrons  versus  their  energy  computed 
for  the  full  absorbed  power  Ws=  10  kW  and  for  different  values  of  Tef 

Figure  1.3.  Upgoing  flux  of  suprathermal  electrons  versus  their  energy.  Solid 
curves  1,  2  and  3  correspond  to  Tef=  5,  7.5  and  10  eV  respectively,  while  the  points  with 
bars  correspond  to  the  observations  made  on  052072. 

Figure  1.4.  Plasma  line  intensities  versus  altitude.  Solid  curves  1,2  and  3 
correspond  to  Tef=  5,  7.5  and  10  eV  respectively,  while  the  points  with  bars  correspond  to 
the  observations  made  on  071392. 
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2.  Optic  emission  from  modified  by  powerful  radio  waves  ionosphere 
Introduction 

Enhancement  of  optic  emission  was  established  in  the  first  ionosphere  modification 
experiments  (Sipler  and  Biondi,  1971;  Adeishvili  et  al.,  1976).  It  was  studied  afterwards  in 
multiple  observations  (Haslett  and  Megill,  1974;  Carlson,  1974;  Bernhardt  et  al.,  1989). 
During  the  last  years  new  interesting  results  have  been  obtained  at  HAARP  facility 
(Peterson,  1999;  Carlson,  1999),  in  sporadic  E-layer  at  Arecibo  and  at  low  duty  circle 
experiments  at  SURA  (Nazyrov  et  al.,  1999).  The  optic  emission  measurements  indicate 
directly  the  existence  of  suprathermal  electrons  accelerated  up  to  the  energies  10  eV  and 
more  in  modified  ionosphere. 

Other  type  of  observations,  demonstrating  the  accelerated  up  to  the  energies  20  eV 
electrons  was  performed  by  Carlson  et  al.,  1982,  Juth  et  al.,  1996,  using  incoherent 
scattering  radar  at  Arecibo.  The  suprathermal  electrons  were  seen  in  these  experiments  in  a 
wide  region  of  the  order  100-km  around  the  reflection  point  of  the  heater  wave.  Some 
indications  of  the  significant  effect  of  suprathermal  electrons  on  SEE  emission  exist  as 
well  (Frolov,  1998). 

Theoretical  considerations  from  the  very  beginning  connected  the  enhancement  of 
optic  emission  with  the  acceleration  of  electrons  in  the  Langmuir  plasma  turbulence  layer 
near  the  reflection  point  of  powerful  HF  O-wave.  Though,  it  was  indicated  as  well,  that  the 
mostly  intensive  red-line  emission  (630  nm)  partly  could  be  determined  by  the  direct 
effective  heating  of  the  main  bulk  of  electrons  (Carlson,  Mantas,  1 996),  or  it's  significant 
part,  captured  inside  nonlinear  density  depletions  -  "striations"  generated  by  the  upper- 
hybrid  turbulence  (Gurevich  and  Milikh,  1998).  The  theory  of  Langmuir  plasma 
turbulence  generated  near  the  reflection  region  of  pump  O-wave  connects  the  electron 
acceleration  process  with  the  modulation  instability  and  with  the  creation  of  nonlinear 
Langmuir  cavitons  -  density  depletions  filled  with  Lagmuir  plasma  oscillations  (Wong  and 
Stenzel,  1975).  The  fast  electrons,  passing  trough  the  turbulent  region  filled  with  cavitons 
obtain  the  energy  from  the  trapped  caviton  plasma  oscillations  in  a  diffusive  way 
(Eastebrook  et  al.,  1975;  Morales  and  Lee,  1977;  Andreev  et  al.,  1980;  Wang,  Goldman 
and  Newman,  1997).  The  main  feature  of  this  mechanism  is  that  the  energy  is  gained  only 
by  fast  electrons  whose  velocity  is  high  enough  v  >  a/co,  where  a  is  a  scale  of  a  caviton,  co 
is  the  heater  wave  plasma  frequency.  Low  energy  electrons  oscillate  in  the  caviton 
adiabatically  and  do  not  get  any  additional  energy.  It  should  be  emphasized,  that  the 
collapse  of  cavitons  does  not  effect  strongly  the  acceleration  process  of  fast  electrons 
(Wang  et  al.,  1997).  In  the  same  time  collapse  can  input  a  significant  part  of  wave  energy 
to  the  thermal  component  (Du  Bois  et  al.  1993,  Hanssen  et  al.  1997).  To  calculate  the  flux 
of  energy,  which  goes  to  fast  electrons  the  detailed,  kinetic  theory  is  needed,  which  takes 
into  account  electron  collisions. 

For  example,  the  number  and  energy  distribution  of  fast  electrons  in  ionosphere 
conditions  is  affected  by  the  so-called  "multiple  acceleration"  (Gurevich  et  al.,  1985, 
Vaskov  et  al.,  1983).  The  multiple  acceleration  process  is  determined  by  the  fact,  that  the 
Langmuir  acceleration  take  place  in  a  strongly  excited  layer,  embedded  in  a  weakly 
disturbed  ionosphere  plasma.  After  the  fast  electrons  leave  the  acceleration  layer  it  collides 
with  a  neutral  molecules  in  nondistrubed  plasma.  Due  to  collisions  some  electrons  could 
return  to  the  acceleration  layer  and  get  there  the  energy  once  more.  The  process  could  be 
repeated  many  times,  thus  creating  a  "multiple  acceleration".  Directly  this  process 
determines  the  high-energy  tail  of  electron  distribution  function  and  a  wide  spread  of  fast 
electrons  around  the  acceleration  region. 
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Though,  understanding  of  the  main  mechanisms  of  the  acceleration  process  in 
Langmuir  turbulence  seems  significantly  clear  enough,  both  the  theory  and  observations  in 
ionosphere  modification  experiments  up  to  now  were  not  developed  to  the  state  when  a 
detailed  quantitative  comparison  between  them  could  take  place.  That  means,  that  some 
fundamental  problems  of  modification  experiments  remain  non-answered.  One  of  them  is 
the  intensity  and  space  distribution  of  the  observed  optic  emission. 

In  the  present  work  we  intend  to  develop  the  existing  theory  and  apply  it  to  obtain 
the  detailed  comparison  with  the  main  part  of  the  existing  experimental  data.  The  paper  is 
constructed  as  follows.  In  section  2  heating  of  electrons  by  both  in  upper  hybrid  and  in 
Langmiur  turbulence  is  considered.  It  is  demonstrated  that  the  calculations  of  electron 
acceleration  in  Langmiur  resonance  layer  without  electron  collisions  lead  to  controversial 
result.  In  section  3  the  theory  of  multiple  acceleration  is  generalized.  The  energy  losses  of 
electrons  on  the  excitation  of  molecular  vibration  level  are  taken  into  account.  In  section  4 
the  model  calculations  describing  the  general  peculiarities  of  the  suprathermal  electron 
space  energy  distribution  are  performed.  The  section  5  describes  the  detailed  comparison 
of  the  theory  with  the  main  observational  data  of  the  artificial  optic  emission  in  F-layer. 
The  strong  enhancement  of  the  optic  emission  in  sporadic  E-layer  modification  is 
discussed  in  sections  6.  In  conclusion  we  formulate  the  nowadays  state  of  the  problem. 
The  strongly  inhomogeneous  character  of  the  disturbed  region  is  stressed. 


1 .2  Heating  of  electrons  in  the  resonance  layer 

In  the  vicinity  of  reflection  point  of  powerful  O-wave  in  ionosphere  both  heating 
and  acceleration  of  plasma  electrons  takes  place. 

As  is  well  known  for  vertically  propagating  pump  O-wave  the  resonance  region  is 
situated  at  the  heights  Z  between  upper  hybrid  Z^  and  Langmuir  Zl  resonances 


^ uh -Z  <  ZL 


0) 


In  this  region  the  frequency  of  the  pump  wave  co  coincide  with  the  frequency  of  the  natural 
oscillations  in  ionosphere  plasma 


4 ne2N{Z, ) 
a>  =  coL  =, - ^ 
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Here  N(Zl)  -  electron  density  at  the  Langmuir  resonance,  where  the  Langmiur  turbulence 
take  place.  And  ZUh  is  defined  by  resonance  condition  for  upper  hybrid  resonance 

ylo>l+0)2c  (3) 
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In  the  vicinity  of  Zuh  the  upper  hybrid  (UH)  turbulence  is  effectively  developing.  The 
width  of  resonance  region  AZ  depends  on  the  pump  wave  frequency  and  plasma  density 
gradient  L: 


AZ  =  Z,  -Z  h  =L-^-r 
L  uh  2o>2 


(4) 


coc  is  electron  gyroferequency.  For  the  ionosphere  F-layer  conditions  usually  AZ  ~  2  +  5 
km. 

Conventional  collision  absorption  of  radio  waves  in  F-region  is  very  weak  -  less 
than  1  dB  [Ginzburg,  1967],  If  plasma  turbulence  in  the  resonance  region  is  excited  the 
absorption  become  abnormally  strong  -  up  to  1CM-20  dB  [Gurevich,  1978]. 

Note,  that  the  Langmuir  turbulent  layer  (LT)  is  placed  in  the  upper  part  of 
resonance  region  close  to  the  reflection  point  (2),  while  close  to  the  upper  hybrid 
resonance  (3)  lays  upper  hybrid  turbulent  layer  (LH).  Acceleration  of  electrons  take  place 
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in  Langmuir  resonance  layer  (we  do  not  speak  here  about  singled  out  case  of  multiple 
gyroresonance  frequency  co  =  ncoc).  The  UH  turbulence  lead  to  formation  of  striations.  An 
effective  anomalous  absorption  of  the  pump  wave  take  place  in  this  region,  which  lies 
below  the  Langmuir  turbulence  layer.  It  means  that  the  UH  anomalous  absorption 
(absorption  on  striations)  can  affect  significantly  the  heating  and  acceleration  processes  in 
the  Langmuir  turbulence  layer  while  diminishing  the  pump-wave  power. 


Heating  of  electrons  in  UH  turbulence  region 

Anomalous  absorption  is  determined  by  the  excitation  of  UH  waves  due  to  direct 
transformation  of  the  pump  wave  on  the  density  gradient  in  striations.  UH  waves  trapped 
in  striations  heat  electrons  and  thus  determine  the  depth  form  and  number  density  of 
striations.  Nonlinear  theory  of  the  electron  heating  and  anomalous  absorption  on  striations 
was  developed  by  Gurevich  et  al.  1996.  The  results  of  the  theory  are  in  a  rather  good 
agreement  with  observations. 

Following  to  the  theory  of  UH  anomalous  absorption  we  determined  the  simple 
relation  for  the  losses  of  pump  wave  power  P  in  the  UH  turbulence  region  on  its  way  to 
Langmuir  turbulent  layer. 


P(ZL)*Po\0 


-K/ 

Xu 
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5.5  +  8,  O<0<1 

Here  K  is  an  attenuation  coefficient  in  dB.  Coefficient  Q  is  proportional  to  the 
relative  depth  of  striations.  It  depends  on  pump  wave  power,  frequency,  plasma 
parameters.  The  state  of  the  striation  development  is  also  significant.  For  the  well 
developed  striations  Q  ~  1 . 


1 .3  Heating  of  electrons  in  Langmuir  turbulent  layer 

The  excitation  of  Langmuir  turbulence,  it  steady  state  and  caviton  formation  in  the 
field  of  powerful  radio  wave  in  conditions  close  to  ionosphere  was  studied  numerically. 
We  will  discuss  the  result  of  the  calculations  of  Hanssen  et  al  and  Wang  et  al  in  order  to 
determine  the  thermal  heating  of  electron  by  Langmuir  turbulence.  The  parameters  used  in 
calculations  are  presented  in  Table  1 . 


Hanssen  et  al 
(1992) 

Wang  et  al  (1997) 

Pump  wave  frequency  (MHz) 

f=  4.01 

f=4.9 

Pump  wave  field  Erms  (V/m) 

1.26 

1.6 

Plasma  density  no  (cm' ) 

ESBflHMHH 

Electron  temperature  (eV) 

0.1 

Pump  wave  energy  Wrms=Erms2/4mioTo 

mmsmumm 

Hanssen  et  al.  (1992)  considered  one-dimensional  Zakharov  system  of  equations 
describing  nonlinear  development  of  plasma  density  and  electric  field  fluctuations  under 
the  action  of  homogeneously  oscillating  electric  field  of  the  pump  wave.  Electron  and  ion 
collisions  and  Landau  damping  were  taken  phenomenologycally,  using  terms  from  linear 
theory.  The  established  average  energy  of  longitudinal  plasma  waves,  concentrated  mostly 
in  cavitons  grew  up  70  times  in  comparison  with  the  pump  wave  energy.  It  gives 
0.19. 

The  power  density  dissipating  by  plasma  oscillations  to  the  heating  of  electrons  is 
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dW  e2E 1  v  n. 

_  e  l 

ds  mco1 

Where  s  is  coordinate  along  the  Earth  magnetic  field,  ve  -  electron  collision 
frequency,  co  -  plasma  wave  frequency  and  Erms2=  4moTe  W1*^.  Following  Hanssen  et  al. 
let  us  suppose  that  the  whole  width  of  the  turbulent  layer  is  of  the  order  of  pump  wave 
length  near  reflecting  point  L  ~  k  «  200  m.  We  obtain  that  the  full  pump  wave  power 
dissipating  in  Langmuir  cavitons  to  the  heating  of  plasma  electrons: 


(6) 


That  value  is  quite  significant.  Really,  if  we  neglect  anomalous  absorption  in 
striations  below  turbulent  layer,  then  to  create  the  electric  field  Erms  =1.26  V/m  (see  Table 
1)  at  the  main  Eiry  maximum  near  reflection  point,  the  pump  wave  power 
60 (sin  a) 4/3 ( 1  OOkm/L ) 1/2  fiK W/m 2  is  needed.  Here  a  is  the  inclination  angle  of  Earth 
magnetic  field.  We  see  that  at  a  ~  1  about  5  dB  of  pump  wave  is  absorbed  due  to  the 
thermal  heating  of  electrons  in  the  cavitons.  Note  that  the  calculated  pump  wave 
absorption  is  effectively  growing  with  latitude  (a4/3).  So,  more  than  60%  of  absorbed 
pump-wave  energy  in  caviton  turbulence  is  going  directly  to  the  thermal  component  of 
electrons. 


Wang  et  al.  1997  elaborated  one-dimensional  kinetic  model.  Not  only  field 
equation  but  also  Vlasov  kinetic  equation  for  noncolliding  electrons  was  integrated 
simultaneously  in  the  model.  The  Langmuir  waves  instability  was  excited  and  established. 
Energy  level  of  average  plasma  oscillations  grew  up  about  30  times  in  comparison  with 
pump  wave  field  reaching  Wrms  =  0.14.  We  see  that  in  spite  of  quite  different  approach,  the 
numerical  value  of  saturated  energy  of  plasma  oscillations  is  close  enough  in  both 
calculations. 

The  power  density  dissipating  to  the  thermal  component  of  electrons  is  again 
determined  by  expression  (6).  The  width  of  turbulent  layer  was  supposed  to  be  A,m  «  180 
m.  From  (6)  we  again  obtain  the  heating  power  Wt=  40  p.  kW/m  ,  what  means  that  about 
2.5  dB  of  pump  wave  energy  goes  to  the  heating  in  cavitons. 

We  see  that  when  Langmuir  turbulence  is  effectively  excited  the  pump  wave 
energy  losses  to  the  heating  of  electrons  is  quite  significant,  compatible  (though  remaining 
less)  to  the  anomalous  absorption  in  UH  layer. 

Significant  advantage  of  Wang  et  al.  1997  kinetic  model  is  the  possibility  to  see 
directly  not  only  saturation  of  plasma  oscillations,  but  acceleration  of  suprathermal 
electrons  as  well.  The  result  of  calculations  shows  that  a  tail  of  initial  Maxwelian 


distribution  function  is  growing  in  time.  An  established  suprathermal  stationary  tail  is 


fs  =8.3xl(T6 


«o. 

VTe 


4<  — <25 

VTe 


(7) 


shown  to  be  approximated  by  power  -  law: 

Here  VTe=(Te/m//2  _  thermal  velocity.  The  one  dimensional  maxwellian  part  of 
distribution  function  has  standard  view: 


fm  = 


»0 


(8) 


The  full  number  of  suprathermal  electrons  is  not  large  about  2%  of  the  thermal  electrons. 
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Let  us  determine  now  the  transport  of  energy  by  fast  electrons.  Integrating 
distribution  function  (7)  we  can  calculate  the  full  power  transported  out  of  the  turbulent 
layer  by  accelerated  suprathermal  electrons 


25  rre  1  ukW 

Wa=  2  J  - mv2vfs dv  *  2.68n0TeVTe  « 1700^- 

J,2  m 


(9) 


We  see  that  the  power  transported  by  suprathermal  electrons  is  20  times  larger, 
than  the  full  power  of  the  pump  wave  (Wp  &  96/jkW/m2).  This  result  is  obviously 
controversial.  It  means,  that  in  spite  of  the  fact  that  mean  free  path  of  the  suprathermal 
electrons  in  ionosphere  F-layer  is  larger  than  the  width  of  the  turbulent  region  Am  *  1 80  m, 
the  electron  collisions  could  not  be  neglected.  The  flux  of  the  energy  out  of  the  layer  is 
determined  by  collisions  _  suprathermal  electrons  not  only  leave  the  layer,  but  due  to 
collisions  they  can  return  back.  Thus  only  the  full  balance  between  accelerated  inside  the 
layer  and  returning  backward  suprathermal  electrons  can  determine  correct  relation 
between  the  pump  wave  energy  losses  and  acceleration.  Just  this  balance  is  studied  in  the 
theory  of  multiple  acceleration.  This  theory  in  modernized  form  would  be  considered 
below. 


Note,  that  the  controversial  result  of  electron  acceleration  in  non  collisional 
kinetics  shows,  that  the  theory  of  plasma  turbulence  saturation  should  be  revised  as  well: 
electron  collisions  outside  the  Langmuir  turbulent  layer  should  be  taken  into  account  to 
organize  a  correct  flow  of  plasma  energy.  In  other  words  boundary  conditions,  which 
connect  turbulent  and  nonturbulent  regions,  could  affect  the  process  of  Langmuir  waves 
saturation  inside  the  turbulent  layer. 


1 .4  Kinetic  theory  of  multiple  acceleration 

To  find  the  distribution  function  we  have  to  solve  kinetic  equation  inside  and 
outside  acceleration  layer  and  match  solutions. 

Inside  the  layer  kinetic  equation  has  a  form: 


df  df  e  „df  _ 

—  +  v„— - E—-0 

dt  dz  m  dv. 


(10) 


Here  E  is  electric  field  of  Langmuir  cavitons  excited  by  powerful  radio  wave 
E(z,t)  =  ^E(z){e-i0*  +c.c) 


The  solution  of  the  equation  (10)  has  the  form: 
f(z,ept)  =  F(z-Z(t),ell  -s(t)) 


(11) 

(12) 


Here  e(t),  Z(t)  the  trajectory  of  a  fast  particle  in  electric  field  E: 


dZ  dV  e  Tr  r— — 

—  =  V,  —  = — E,  V  =  f2  elm 

dt  dt  m 

To  solve  our  problem  we  need  to  determine  the  changing  of  distribution  function 
Af  when  a  fast  particle  crossing  the  acceleration  layer.  Averaging  relation  (12)  on  fast 
oscillations  we  find  the  changing  of  distribution  function  after  passing  through  the 
acceleration  layer). 


v-i-l 

2  ds„ 


SF_ 

de„ 


(13) 


According  to  (11)  after  averaging  on  fast  oscillations  in  a  fix  spatial  point  energy  gain  is 
given  by  a  formulae: 
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h.s{z,i)  =  -Afcos 


(14) 


+00 

As(^t)=^Ek\,  Ek=  jE(z)e~,kzdz, 


For  example,  for  Gauss  like  caviton  E(z)=  Eoexp(-(z/a)2)  we  have  from  (11): 


A.e  =  n'2eE^a  exp(-  ka 2 ) 


(15) 


We  see  from  (15)  that  only  electrons  with  high  enough  energy 


=  m(coa)2  =  Te 


(  \ 
a 

2 

fa] 

UJ 

l»e) 

»1 


(16) 


where  De  is  a  Debay  length,  could  be  accelerated.  It  is  easy  to  see  that  relation  (16)  has  a 
universal  character  and  does  not  depends  on  a  real  form  of  a  caviton.  It  is  obvious  that  if 
we  have  some  number  of  noncorrelated  cavitons  the  average  energy  changing  is  equal  to  a 
sum 

(A*2)=Z(A  £i) 

i 

Now  let  us  consider  the  solution  outside  acceleration  layer  at  z>0  and  z<0.  Kinetic 
equation  has  a  form: 


df  df  c 
—  +  vu—  =  ~S 
dt  dz 


(17) 


Here  f(t,v,fj,z )  is  the  distribution  function  of  electrons,  S  -  collision  integral,  p  _ 
cosine  of  an  angle  between  the  velocity  of  electron  and  z  direction.  The  collision  integral 
consists  of  two  parts:  inelastic  So  and  elastic  Si.  The  collision  integral  depends  strongly  on 
energy.  For  energies  e<s*  (s*«4  eV)  the  main  collision  crossection  is  defined  by 
oscillational  levels  and  inelastic  collision  integral  has  a  diffusive  form.  But  for  energies 
e>£*  the  collision  crossection  is  determined  by  optical  levels.  Thus  the  collision  integral 
could  be  presented  in  a  form: 

'  df^Qf' 


S0(F0)  =  - 


de 


R- 


de 


S0(F0)  =  v0(e)F0 


e  <  e 


V0  =  V°0Nn, 


S  >  £ 


(18) 


V0=VCT,Wm 


5,(/1)  =  v1(f)/i, 

Srt=<T,0(e-et  )+%&(£-£*) 

t 

The  oo  is  the  sum  of  effective  crossection  of  inelastic  collisions  with  neutral  atoms  Nm  -  the 
total  density  of  neutrals.  at  is  transport  electron  crossection  of  electron  elastic  collisions 


O’o 


N, 


mk 


N„. 


(19) 


The  averaged  part  of  total  energy  of  an  electron  loosing  in  one  collision  is  small;  it  is 
defined  by  parameter 


S  =  vjv{  =ajat 


(20) 
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High  electron  energies(e  >  £*) 

Because  of  (20)  we  can  expand  the  equation  (17)  on  small  parameter  5.  At  first 
approximation  the  distribution  function  does  not  depend  on  p.  Taking  into  account  first 
approximation  one  can  find: 


df} 

dt 


o  +y¥i 


3  dz 


■-sQ, 

dFg 

dz 


(21) 


:  =  -vJi 


Boring  in  mined  the  boundary  conditions  (13)  at  z  =  0  in  stationary  conditions  from  (21) 
one  can  find: 


9j -  =  sign(z)^^-\  J(Ae2 (ff, ))<fe, 


dz 


9vs  ds 


(22) 


It  is  easy  to  see  that  for  collision  integrals  in  the  form  (18)  the  equation  (22)  depends  on 
energy  parametrically  only.  So  we  can  rewrite  equation  (1 8)  in  the  region  of  energies 

^r  =  3  SF0,  $=[^dz=[o{Nm(z)dz  (23) 


Taking  into  account  that  8  is  a  slow  changing  function  it  is  possible  to  solve  (23)  in  WKB 
approximation 


F0(£,z)  =  F0(e) 


(24) 


/i  (£ ,  z)  =  sign(z)(3S )^  F0  (£,  z) 

Here  Le  is  a  relaxation  scale  of  fast  electrons  with  given  energy 


1 

Wm(30-,cro)>2 


Middle  electron  energies  (Te<£<£*) 

In  the  region  of  energies  £<  £*  the  kinetic  equation  (21)  takes  the  form: 


d2Fn 


d d£ 


sIrSF, 
d£  j 


(25) 


Since  the  solution  of  kinetic  equation  in  the  region  £>  e*  has  a  form  (24),  to  match  this 
solution  with  solution  in  region  £<  £*  we  had  to  search  the  solution  for  e<  £*  in  a  form 


f 

zr  dz 

1 

lk{z)\ 

F0  =  g(f)exp|- 

The  equation  (25)  takes  the  form: 

de  \  d£ )  vj 
Taking  into  account  that 

v0v,  «  Rg 

One  can  see  from  (26)  that  at  first  approximation 
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£<£ 


F0(e)  =  F0(£*)  =  const  £<£  (27) 

The  distribution  function  Fo(s)=F0(s,z)|z=o  according  to  boundary  condition  (13)  is 
determined  by  diffusion  equation: 


and  F0(s)  =  F0(e’)  =  const,  e<e * 

We  introduce  here  effective  temperature : 

T€  (e)  =  kKe,  A  e  =  j-  j(Ae2  (fj  ))<fc, 


F0=f,-f(0\s) 


X 


3X  _y 

K  =  —i=8  A 


2V2 


(28) 


(29) 


here  As  is  mean  electron  energy  increasing  in  accelerating  layer.  The  effective  temperature 
increase  gradually  with  increasing  energy  and  saturate  at  the  level. 

Teffm  =K(n/2/2&£m,  sm<£, 

Here  n  is  the  number  of  cavitons  in  accelerating  layer. 

The  intensity  of  electron  acceleration  is  defined  by  parameter 


eff 


T  =  -fm(£) 


dfm(£ ) 


T  '  '  '  7  ds 

For  weak  acceleration  (y<l)  the  distribution  function  of  electrons  does  not  change 
substantially 


F0(e)  = 


\  T  , 


/<”>(«) 


In  opposite  case  of  strong  acceleration  y  »  1  the  equation  (29)  in  WKB  approximation 
has  the  following  solution: 

fo(£)  =  fm(£)  £<F 


Ih 


Me)  =  foM 


X ,  .  _i/ 

M  exp 

1 - 

1 _ ! 

2  £Teffm(£)Sy\£) 

1  *  J 

3 

«o" 

II 

_ 1 

£>£„ 


(30) 


The  value  sth  in  (30)  is  defined  by  matching  of  distribution  functions  in  saddle  point 
d'effm  ~  T(£lh  ). 

;l/2 


In  the  case  of  strong  acceleration  in  the  limit  Teff  ~  const,  5  «  const  the  distribution 

function  (30)  has  the  form: 


m=cK0 


r  E^ 


(31) 


Here  Ko(x)  is  a  modified  Bessel  function.  The  number  of  fast  electrons  is  defined  by  the 
integral 

N{  =  jf(s)d\ 


But  the  normalized  constant  C  is  determined  by  the  microwave  power  absorbed  by  fast 
electrons.  The  absorbed  power  consists  of  two  terms 
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(32) 


P  =  Pl+P2 


P  - i6n 

2~  m2  VI 


C\e> 


£,h 


\T<#  J 


The  first  term  Pi  in  (32)  is  responsible  for  loses  on  rotational  levels  of  molecules  and  the 
second  one  P2  describes  the  electron  acceleration. 
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4.  Drift  effects  on  the  scale  of  field-aligned  irregularities 


Irregularities  of  electron  temperature  and  plasma  density  in  the  F-region  of  the  ionosphere 
are  strongly  elongated  magnetic  field  due  to  dramatic  difference  in  transport  coefficients 
along  the  magnetic  field  and  transverse  to  it.  Their  scales  along  the  magnetic  field  and 
reply  to  heating  sources  in  stationary  state  were  determined  previously  neglecting  plasma 
drift  effects  both  in  linear  [Gurevich,  1978]  and  many  other  publications]  and  nonlinear 
[Gurevich  et  al.,  1995]  approximations.  These  scales  are  rather  long  reaching  several  ten 
kilometers.  Everybody  settle  that  such  elongated  irregularities  should  be  affected  by 
ionospheric  drifts  perpendicular  to  magnetic  field,  but  up  to  date  there  was  no  theoretical 
estimations  for  this  influence.  It  can  be  just  supposed  that  drifts  could  substantially 
diminish  field-aligned  scales  and  efficiency  of  excitation  of  irregularities  with  sufficiently 
small  scales  in  the  drift  direction.  Here  we  will  attempt  to  consider  drift  influence  on 
electron  temperature  and  plasma  density  irregularities  under  approximation  that  source 
dimension  along  magnetic  field  is  much  less  than  characteristic  transport  lengths,  and 
neglecting  all  transport  processes  across  magnetic  field.  Two  main  problems  will  be 
considered:  constant  drift  effects  on  irregularities  caused  by  fixed  in  the  ionosphere 
sources  (for  example,  narrow  radio  beam),  and  gradient  drift  effects  on  irregularities 
created  by  moving  with  plasma  sources  (like  excitation  of  irregularities  by  plain 
electromagnetic  wave). 


Constant  drift  velocity 

Transport  equations  along  the  magnetic  field  for  disturbances  of  electron 
temperature  T  and  number  density  N  can  be  written  taking  into  account  homogeneous  drift 
with  constant  velocity  Vo  across  magnetic  field  in  the  form  [Gurevich,  1978] 


-v>j--*r  =  -QXz) 

OZ  OX 


(1) 


A 


d2N 

dz2 


dN_ 

dx 


( 


-yN  =  - 


QN 


D 


m 


K„ 


5{z)- 


D. 


T\  I 


ar 

dx 


+  8vT 


Here  D\\,  Dj]\,  and  kjj  are  diffusion,  thermal  diffusion,  and  heat  conductivity  along  the 
magnetic  field  respectively,  y  is  the  recombination  coefficient,  8  is  the  average  part  of 
Selection  energy  loss  under  collisions  with  ions  and  neutrals,  and  v  is  the  election  collision 
frequency.  It  is  supposed  in  (1)  that  irregularities  are  induced  by  fixed  heating  Qt  and 
density  Qn  sources  which  scales  along  the  magnetic  field  is  much  less  than  characteristic 
scales  of  transport  processes  so  that  Dirac  delta  function  8{z)  can  be  used.  Coordinate  z  is 
along  the  magnetic  field,  and  x  is  along  drift  velocity  Vo  perpendicular  to  the  magnetic 
field. 

Introducing  un-dimensional  variables  according  to 

T  N 

t  =  — ,  n  =  ■ — - 

Tn  N, 


lo  o 

with  characteristic  lengths 


£  =  — 
’  L  ’ 

X,y 


Z 

^  11' 


(2) 


L  =Ll 

v  Sv 

and  un-dimensional  heating  and  density  sources  as 


(3) 
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qT  = 


1 


9 "  1 V„^C, 

where  To  and  A'o  are  undisturbed  electron  temperature  and  number  density  respectively, 
and  looking  for  a  solution  of  (1)  in  the  form  (x,n)  cc  exp(ik£),  one  can  rewrite  transport 
equations  (1)  as 

d  z 


'Qn> 


(4) 


~(/£  +  l)r  =  £(^), 

-  (/A  +  p2  )n  =  ~(qN  -  RqT  )S  (g )  -  R(ik  +  l)r 


d£ 

d2n 

dg 2 

Following  denotations  are  used  in  (5): 


T  D 

ft  _  1o1jt\\ 


No  A 


(6) 


Solution  of  (5)  depends  on  |£|  (that  can  be  proved  by  substitution  £->  -Q,  and  only  a 
region  of  £>  0  will  be  considered.  In  this  case  (5)  takes  a  form  of 

~-(ik  +  l)r  =  0,  (7) 

d£ 

d2n 

dg 2 

with  boundary  conditions 
dr 


-  ( ik  +  l)r  =  0, 

-  {ij3k  +  fj2)n  =  -R(ik  +  l)r 


I 


n 


Ig*— >0O 


=  0’  2^Uo  =-qr, 

=  0,  2^\  0=-(qN-RqT). 

d£ 


(8) 


Solution  of  (7),  (8)  is  straightforward  and  leads  to 

1  -«  1  -if  1 

t  =  —  qTe  ,  n  =  — qTe  4 
2  />  2  s 


Rqr{p' 


>e-p<  -se'5* 


\ 


where  p  =  -J\+~ik  and  s  =  y]p2  +  ifik  and  Kip),  7?(  s)  >  0  H  stands  for  the  real  part). 

Let  us  analyze  (9)  supposing  for  simplicity  q n  =  0.  Amplitudes  of  electron  temperature  and 
density  irregularities  produced  by  constant  heating  source  qx  depend  on  its  scale  along 
drift  direction  decreasing  for  smaller  scales.  The  same  dependence  exists  for  the  length  of 
irregularities  along  the  magnetic  field.  There  is  also  phase  shift  between  heating  source 
and  irregularities.  Temperature  disturbances  can  be  characterized  by  "efficiency"  of 
heating  Gv  and  their  relative  to  heating  sources  phases  @T  that  can  be  defined  as 

l  ,,  _  ,  . 

=  arg(r0)  =  arg 


(9) 


_2_ 

qT 


Gr  =  —  r0  =  ■  , - 

1  |VT Tiii 


=  0  +  k2)Y\ 


—f==  =  arctan  k.  (10) 
■ql  +  ik  2 


Similar  parameters  of  "efficiency"  G„  and  phase  On  can  be  defined  for  plasma  density 
disturbances 

2  i  i  1  +  p  ,  1 


G„  =■ 


=  arg(«0)  =  arg- 


RRt  '  '  I  Jl  +  ik  +  yjp2  +  ifik  I  ”  Vl  +  ik  +  y/ju2  +  ifik 

Characteristic  length  Ax  of  heated  region  along  magnetic  field  can  be  estimated  as 


(11) 
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(12) 


@r  =  ~\T0 


41 

dq 


Jg->  0 


vj 


1  +  Vl  +  yP 


Similarly,  field  aligned  scale  of  plasma  density  irregularities  An  can  be  estimated  from 
f  jl  jv* 


©„  =  -In, 


G 


d  n 

dg 


(13) 


Jg^O 


In  the  long  scale  limit  k->0  one  can  obtain  that 

1,  G=— — ,  <£  =<D  =0, 

"  \  +  m 


@,  =  1 


(14) 


as  in  case  of  absence  of  drifts.  It  is  obvious  that  drift  does  not  affect  disturbances  with  long 
enough  scales  along  drift  velocity. 

In  the  short  scale  limit  k—>°o  there  is  significant  dependence  of  "efficiencies"  and 
characteristic  scales  of  irregularities  along  the  magnetic  field  on  wave  vector  k : 

1  „  1  x  ^  l~2 

g„  =-f — P=T7=,  a>T=<s> 


G  = 


4k’ 


X+JpW' 


a  = 


’’Efficiencies”  of  excitation  of  irregularities  for  short  scales  decrease  as  fcvl  as  well  as  their 
characteristic  lengths  along  magnetic  field.  Phases  of  both  electron  temperature  and 
density  disturbances  become  constant  in  quadrate  to  heating  source. 

Dependencies  of  above  parameters  on  wave  number  k  are  shown  in  Figure  1.  Contours  of 
constant  temperature  and  plasma  density  are  shown  in  Figure  2  for  some  values  of  wave 
number  k. 

3.3  Constant  drift  gradient 

Transport  equations  along  the  magnetic  field  for  disturbances  of  electron  temperature  T 
and  number  density  N  in  the  coordinate  frame  connected  to  moving  with  plasma  sources 
can  be  written  taking  into  account  constant  drift  gradient  dV(/dx{z}  =const  in  the  form 
[Gurevich,  1978] 


KTi, 


A 


d2T  dVn  dT 


dz2 
d2N 
dz 2 


-z- 


dz  dx 


-SvT  =  -QrS(z\ 


dV0  dN 
dz  dx 


■yN  =  - 


A. 


Qn—^Q, 


Dt 


S(z)-2L 1 


II  V 


dV0  ST 
dz  dx 


+  SvT 


(16) 


J 


Used  notation  is  the  same  as  in  (1).  It  is  supposed  that  irregularities  are  induced  by  heating 
Qt  and  density  Qn  sources  that  move  with  plasma  drift  velocity  Vo. 

Un-dimensional  equations  for  variables  (2)  where  Ly  (3)  is  substituted  by 


L  =  Lt  dV0 
v  Sv  dz 


(17) 


for  harmonic  solutions  of  (16)  (z,n)oc  exp(ikd)  can  be  written  as 


^-(ik£  +  \)v  =  -qT8{g\ 

dg 

~j  -  {ifiC  +  m2)"  =  ~{qN  -  RqMs)  -  R(ikC  + 

dq 

Un-dimensional  sources  (4)  and  denotations  (6)  are  used  here. 


(18) 


(15) 
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Substitution  £  -»  to  (18)  leads  to  complex  conjugate  system  of  equations  that  means  that 
(z,n)  (-Q  =  (r,n)  (Q.  Considering  region  of  g>  0  one  can  rewrite  (18)  as 
d2r 


dC 

d2n 


2  (ik£  +  l)r  =  0, 

-  (i/3k£  +  /j2)n  =  -R(ik£  +  l)r 


(19) 


dg2 

with  boundary  conditions 
^  dr 


rU=o. 


29? 


\g->0 


"U«=0: 


29? 


dt; 

dn  . 
~d~^ 


Qt  >  ^  7  I (=0~  ® 


(20) 


i  <r=o‘ 


0 


lN~Rclr\  3«| 

J 

where  9?  and  3  stand  for  real  and  imaginary  parts  respectively.  Relations  for  imaginary 
parts  of  x  and  n  at  {=  0  are  followed  from  continuity  conditions  at  the  source. 

Solution  of  equation  (19)  for  x  with  appropriate  boundary  conditions  (20)  is 


7  =  -^-Ai,[v(0)]Ai[v(4')], 


v(<r) =(/*)%+ T 


ik 


(21) 


where  Ai(z)  is  Airy  function  that  does  not  diverge  at  infinity,  and 
Ar  =  9?{Ai[v(0)]}9?|/^)^  Ai  [v(0)]}+  3{Ai[v(o)]}3^)^  Ai'  [v(o)]}  (22) 

Substituting  (21)  into  density  equation  (19)  and  solving  it  with  appropriated  boundary 
conditions  (20)  one  can  obtain 

»  =  -Tf-Ai,[w(°)]Ai[*fe')]+T?LArM0)]Ai[»(f)] 


2  A., 


2A„ 


■3 


3 


1  +  i«|tt)^JAi-[v(0)]Bi'[v(0)]} 

Ai'[w(0)f  AiMf)] 
J  Ai*  [v(o)]B  i[w(o)] 

B  i  [w  (0  )]  v  (^  )  A  i  [v  (^  )]  A  i  [w  (^  )] 

( 

Ai[w(4-)]Jrff'v(f)Ai[v(f)]Bi[w(f )] 


2A 
ik 

Jj 


( ik  \'A 

f 


(23) 


w  =  (/#H  £  +  -£- 
[  ifik) 

where  Bi (z)  is  another  Airy  function.  Here  in  (23) 

A„  =  9?{Ai [w(o)]}9? \}PkY3  Ai’[w(o)]}+  3{Ai[w(o)]}3f//ft)^  Ai'[w(o)]} 
and 


(24) 
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J  =  ]^'v(UMvfc  )Mwt’)l 
0 

Using  asymptotic  expansions  of  Airy  functions  for  \z\  —>co 

Ai(z)  ->  *  y  exp  j-  \z/2\ 

2-v/^r z'4  l  3  J 


Bi(z)  -» 


2 


I2  3A\  ■ 

exp  j~z  ^ -/exp 


one  can  show  that  in  the  long  scale  limit  obtained  expressions  (21)  and  (23)  for  rand  n 
respectively  can  be  reduced  to  well  known  expressions  for  the  case  without  drifts: 


2  Rq  ' 


Tk->0  =  \qTe<>  «*-*>  =\q»e~K  Rq‘  ^  ~ *** )  {21) 

Suppose  now  for  simplicity  that  there  is  no  source  of  plasma  density  qu  =  0.  In  this  case 
phases  of  temperature  and  density  disturbances  coincide  with  the  phase  of  heating  source 
due  to  relation  (n,  t)(-Q  =  (n,  r)*( Q. 

“Efficiencies”'  of  heating  Gz  and  excitation  of  density  irregularities  G„  are  defined  by  (10) 
and  (11),  while  characteristic  lengths  along  magnetic  field  of  temperature  Lz  and  density 
L,,  irregularities  —  by  (12)  and  (13)  respectively.  It  can  be  shown  that  GT=  Lr. 

It  is  obvious  that  in  the  long  scale  limit  £-»0. 

o, = 4 = i.  g"=Tw;’  (28) 

In  this  case  drift  does  not  affect  development  of  irregularities. 

In  the  short  scale  limit  k~> oo 

,29) 

V3  Ai  (0)  V3  Ai  (0)  i  +  pAi  +  pA 

It  is  taken  into  account  here  that 


Ai(0)Ai'(0)J 


[dw-MAi(w)Ai'(  =  -■  1  ~ - y- 

o  V  J  \  +  pyi+pyi 


Let  us  also  note  that 


2  Ai(o)  _  2r\}j. 

V3  Ai  (o)  3%  y{2/ 


-1.58393 


where  IJx)  is  Gamma  function.  "Efficiencies"  of  excitation  of  irregularities  for  short 
scales  decrease  as  K>/1  as  well  as  their  characteristics  lengths  along  magnetic  field. 
Dependencies  of ' '  efficiencies"  Gx  and  Gn  on  wave  number  k  are  shown  in  Figure  3 . 


3.4  Conclusions 

Let  us  estimate  relevance  of  above  considerations  to  the  Earth's  ionosphere.  Typical  values 
of  ionospheric  parameters  in  the  F-layer  are  (Sv)'1  =  10-30  s,  Lt  =30-100  km,  and  drift 
velocity  at  mid  latitudes  can  vary  in  the  range  of  Vo  =  10-200  m/s.  As  it  can  be  seen  from 
Figures  1,  3  drifts  can  essentially  affect  development  of  irregularities  with  scales  of  order 
or  less  than  Ly  for  fixed  heating  source  and  Ly  for  a  source  moving  with  plasma  drift. 
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In  case  of  ionospheric  modification  by  powerful  radio  waves  fixed  source  is  realized  for 
narrow  beams,  for  radio  waves  focused  by  natural  horizontal  inhomogeneities  in  lower 
layers  of  ionosphere,  and  for  self-focusing  of  powerful  radiation.  Typical  values  of 
characteristic  scales  of  irregularities  when  drift  influence  becomes  significant  are  Ly  = 
100-6000  m  for  parameters  given  above.  These  scales  lay  in  typical  interval  both  for 
natural  focusing  and  for  self-focusing  of  radio  waves. 

Drift  velocity  gradient  becomes  essential  when  assumption  of  homogeneous  electric  field 
of  powerful  radio  wave  is  applicable  for  description  of  excitation  of  irregularities  as  in 
case  of  resonance  instability  or  thermal  parametric  instability.  Gradient  of  drift  velocity 
can  exist  in  the  ionosphere  even  just  due  to  the  dependence  of  magnetic  field  strength  on 
the  altitude.  Drift  velocity  in  crossed  electric  E  and  magnetic  B  fields  is  described  by 
relation  Vo  =  cE/B.  Supposing  that  magnetic  field  lines  are  equi-potential  one  can  found 
that  EocS~!/~  where  S  is  area  of  magnetic  flux  tube,  Sxl/B.  Hence,  Vo  xB'1/2.  Considering 
Earth's  magnetic  field  as  dipole  field  with  B  x  R'3  one  can  obtain  for  drift  velocity  Vo  x 
R3/2  where  R  is  distance  from  the  center  of  the  Earth,  and  dVo/dz  =3  V(/2R.  Gradient  of 
drift  velocity  in  the  F-region  of  the  ionosphere  can  be  estimated  as  0.2-4. 5  sm/s  /km  for 
above  parameters.  Characteristic  scales  for  drift  gradient  Ly  =  0.6-150  m  are  comparable 
to  those  usually  considered  in  theories  of  small  scale  stratification  of  ionospheric  plasma 
by  powerful  radio  waves.  Drift  gradient  can  also  essentially  curve  irregularities  from  the 
straight  line  along  magnetic  field  and  determine  aspect  sensitivity  of  radio  wave  scattering 
off. 

References 

A.V.  Gurevich,  Nonlinear  phenomena  in  the  ionosphere,  Springer,  Berlin,  1978. 

A.V.  Gurevich,  A.V.  Lukyanov,  K.P.  Zybin,  Stationary  state  of  isolated  striations 
developed  during  ionosphere  modification,  Phys.  Lett.  A,  206,  247-259,  1995. 


24 


Figure  3.1:  "Efficiencies"  of  heating  GT  (a)  and  generation  of  irregularities  Gn(d),  phases 
of  heating  <f>T  (b)  and  irregularities  On  (e),  and  characteristic  lengths  along  magnetic  field 
of  electron  temperature  AT(c)  and  density  irregularities  An(f)  versus  transverse  scales  of 
irregularities.  In  panels  (d)  ,  (e),  and  (f)  curves  are  shown  for  J3=l  and  /r=0  (solid  lines), 
0.3  (dots),  1  (dashes),  3  (dash-dots). 
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Figure  3.2:  Contours  of  constant  electron  temperature  /(a),  (b),  (c)/  and  plasma  density 
/(d),  (e),  (f)/  for  fJ=\  and  ju=0  for  different  values  of  wave  number  k.  Along  horizontal  axis: 
4  =  x/Ly,  along  vertical  axis:  4=  z/Lj.  Panels  (a)  and  (d)  correspond  to  kLy  =  0.1,  absolute 
contour  values  for  panel  (a)  are  0.9,  0.6,  0.3,  0,  for  panel  (d)  —  0.7,  0.5,  0.3,  0.  Panels  (b) 
and  (e)  correspond  to  $kL_{V}  =  1$,  absolute  contour  values  for  panel  (b)  are  0.7,  0.5, 
0.3,  0,  for  panel  (e)  —  0.4,  0.3,  0.2,  0.  Panels  (c)  and  (f)  correspond  to  kLv  =  10,  absolute 
contour  values  for  panel  (c)  are  0.25,  0.2,  0.15,  0,  for  panel  (f)  —  0.15,  0.1 1,  0.07,  0. 
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Figure  3.3:  "Efficiencies"  of  heating  Gr  =  AT  (a)  and  generation  of  irregularities  G„  (b) 
versus  transverse  scales  of  irregularities.  Curves  in  panel  (b)  are  shown  for  p  =  1  and  p  = 
0,  0.3,  1,  3  (lower  curve  corresponds  to  n=0). 
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The  work  is  driven  by  the  intersection  of  two  lines  of  activity,  one  experimental, 
the  other  theoretical,  that  converge  on  an  importantly  new  view  of  the  physics  of  how  high 
power  HF  radio  waves  interact  with  collisional  plasmas,  such  as  the  earth's  ionosphere. 

Its  goal  is  to  seize  this  opportunity  to  amalgamate  the  complementary  US  and 
Russian  recent  findings,  to  strengthen  the  foundation  of  plasma  physics  underlying  future 
work  in  the  growing  field  of  ionospheric  HF  heating.  The  following  main  problems  are 
considered  in  the  present  Interim  Report: 

1 .  Observations  of  fast  electrons  in  a  wide  region  around  O-wave  reflection  point 
and  comparison  of  observational  data  with  multiple  acceleration  theory. 

2.  Energization  of  electrons  by  Langmuir  turbulence  in  the  plasma  resonance 
region  and  its  relation  to  the  observations  of  excited  optical  emission. 

3.  Drift  effects  on  the  scale  of  field-aligned  irregularities. 
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1.  Experimental  Studies  of  Energetic  Electrons  in  Ionosphere  Using  ISR 
Technique  and  Comparison  of  Observational  Data  with  Multiple  Acceleration 

Theory 


Abstract 

Two  sets  of  observations  of  suprathermal  electrons,  produced  by  the  interaction  of 
powerful  radio  wave  with  the  ionosphere,  as  seen  using  the  incoherent  scatter  radar  (ISR) 
technique  are  presented.  The  observational  data  are  compared  with  the  theory  of  multiple 
acceleration  of  electrons  in  the  strongly  excited  resonance  region  near  the  reflection  point 
of  the  powerful  radio  wave.  The  structure  of  the  wide  perturbed  region  filled  with 
energetic  (10-20)  eV  electrons  is  determined.  The  size  of  this  region  along  the  Earth 
magnetic  field  is  shown  to  be  1 00  km.  The  full  power  going  to  the  accelerated  particles  is 
determined:  it  is  6-8  kW  or  (4-6)%  of  the  entire  HF  radiated  power.  The  power  carried  by 
the  suprathermal  electrons  flux  escaping  into  magnetosphere  is  of  the  order  1  kW. 

1.1  Introduction 

Observations  of  fast  suprathermal  electrons,  due  to  the  interaction  of  powerful 
radio  waves  with  the  ionosphere,  using  incoherent  scatter  radar  (ISR)  technique  were  first 
proposed  and  realized  by  Carlson,  Wickwar  and  Mantas  in  1972  (see  Carlson  et  al. 
[1982]).  The  powerful  O-mode  wave  in  the  resonance  region  near  the  reflection  point 
effectively  excites  plasma  waves  -  natural  oscillations  of  the  ionosphere  plasma.  Plasma 
waves  due  to  the  nonlinear  cavitation  process  accelerate  suprathermal  electrons  to  energies 
in  excess  off  -  10+20  eV,  i.e.  two  orders  of  magnitude  higher  than  the  thermal  electron 
energy  Tc»  0.1  eV.  The  fast  electrons  propagate  in  the  ionosphere  to  large  distances,  10- 
100  km  from  the  acceleration  region.  Those  electrons  excite  nonthermal  plasma  waves, 
which  were  detected  by  ISR  [Carlson  et  al.,  1982], 

We  discuss  here  the  results  of  two  sets  of  experiments.  The  experimental  data  are 
here  compared  in  detail  with  theory  based  on  the  mechanism  of  multiple  electron 
acceleration. 

Carlson  et  al.  [1982]  gave  a  qualitative  explanation  of  their  observations.  They 
pointed  out  the  significant  role  played  by  multiple  scattering  of  electrons  in  the  neutral 
atmosphere.  The  quantitative  theory  of  multiple  acceleration  was  proposed  by  Gurevich  et 
al.  [1985].  The  comparison  of  the  theory  and  observations  can  help  reach  a  still  higher 
level  in  our  understanding  of  a  number  of  space  plasma  phenomena  in  which  accelerated 
electrons  play  an  important  role. 


1.2  Observations 

The  most  sensitive  and  direct  ground  based  means  of  detecting  10+20  eV 
suprathermal  electrons  is  by  incoherent  scatter  plasma  line  observations  of  the  Langmuir 
waves  produced  in  the  ionosphere  by  these  electrons  [Perkins  and  Salpeter,  1965].  This 
technique  has  been  used  to  study  photoelectron  fluxes  in  the  ionosphere  [Yngvesson  and 
Perkins,  1968].  Carlson  et  al.  [1982]  applied  the  ISR  technique  to  measurements  of 
suprathermal  electrons  produced  by  the  interaction  of  powerful  radio  wave  with  the 
ionosphere  at  Arecibo.  We  present  here  two  types  of  experimental  data.  The  primary  and 
secondary  data  were  gathered  on  the  nights  20  May  1972  and  13  July  1992  respectively. 
The  primary  data  were  published  by  Carlson  et  al.  [1982],  here  we  will  revisit  those  results 
and  supplement  them  with  the  complimentary  set  of  secondary  1992  data  [Djuth  et  al., 
1996], 
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For  the  1972  Arecibo  data  here,  a  HF  dipole  feed  hung  over  the  1000-foot  diameter 
ISR  radar  dishes,  for  which  the  138  kW  transmitter  at  7.63  MHz  delivered  about  30  W/m2 
to  the  HF  reflection  height.  For  the  1992  Arecibo  data  here,  a  separate  48  array  of  log 
periodic  antennas  17  km  away  from  the  ISR  radar  dish,  for  which  the  300  kW  transmitters 
at  5.1  MHz,  into  50%  efficiency  transmission  lines,  delivered  about  40  W/m2  to  the  HF 
reflection  height. 

For  20  May  1972,  strong  plasma  waves  were  excited  near  the  reflection  point  at  the 
altitude  of  z  =  285  km.  Night  time  plasma  line  intensities  were  observed  to  be  enhanced  by 
a  factor  of  10-100,  over  a  range  of  altitudes  which  extended  to  below  250  km.  They 
relaxed  to  their  normal  level  within  the  electron  transport  time  after  the  transmitter  was 
turned  off. 

The  430  MHz  incoherent  scatter  radar  was  used  to  diagnose  the  ionosphere,  which 
allowed  measurement  of  altitude  profiles  of  the  background  plasma  temperature  Te  and 
electron  concentration  N.  The  plasma  line  echo  from  the  plasma  waves  excited  by  the 
suprathermal  electrons  provided  the  altitude  of  a  set  of  plasma  frequencies  fp  between  5 
and  7.5  MHz,  with  an  altitude  resolution  of  1.5  km. 

For  a  given  diagnostic  radar  wavelength  Ar  the  main  scattered  signal  comes  from 
ionosphere  plasma  waves,  whose  wave  vector  is  directed  toward  the  radar,  and  the  phase 
velocity  is  equal  to  the  velocity  of  suprathermal  electrons  v 

v„=Tt,4  0) 

Here  fp  =  (e2N/47rm)l/2  is  the  local  plasma  frequency.  For  the  Arecibo  radar  Ar  =  70 
cm  and  from  eq.  (1)  one  can  obtain  that  the  energy  of  electrons  is 

e  -  ^mv2ph  =  0.35/^F  (2) 

The  intensities  of  plasma  line  echoes  were  shown  in  Fig.l  of  Carlson  et  al.  [1982], 
A  significant  enhancement  of  radar  echo  intensities  is  evident  at  large  distance  below  the 
reflection  point  of  the  radio  wave  at  285  km.  The  upgoing  one-dimensional  fluxes  of 
suprathermal  electrons  at  the  heights  266  km  and  256  km  determined  from  these  data  are 
shown  in  Fig.3  of  Carlson  et  al.  [1982]. 

The  additional  data  were  gathered  on  the  night  13  July  1992.  The  reflection  altitude 
was  295  km,  with  strong  enhancements  of  plasma  line  intensity  observed  in  the  vicinity  of 
the  reflection  level.  An  example  of  plasma  line  echoes  collected  at  the  heights  337-367 
km,  where  the  ISR  beam  intersected  the  flux  of  the  accelerated  electrons,  is  shown  in 
Fig.  1.1.  The  echoes  were  obtained  both  below  and  above  the  maximum  of  F-layer. 
Electron  density  Nmax~  5.5- 1 05  cm'3  was  reached  at  F-maximum  heights  z  =  350  km.  The 
observed  plasma  line  echoes  correspond  to  upgoing  electron  flux.  Note  particularly  here  a 
significant  enhancement  of  plasma  line  echoes  even  at  altitudes  more  than  70  km  above 
the  reflection  layer. 


1 .3  Brief  Outline  of  the  Theory 

In  ionosphere  modification  experiments,  electrons  gain  energy  in  a  strongly 
disturbed  Langmuir  resonance  layer  near  the  reflection  point  of  the  O-mode  wave.  The 
essential  point  of  multiple  acceleration  theory  is  a  large  fraction  of  fast  electrons,  after 
leaving  the  acceleration  layer  will  return  back  due  to  collisions  with  neutral  particles,  and 
thus  can  gain  additional  energy.  The  process  can  be  repeated  many  times.  Thereby,  a  wide 
region  around  the  acceleration  layer,  which  is  elongated  along  the  geomagnetic  field,  can 
be  filled  with  strongly  heated  magnetized  suprathermal  electrons.  The  elongation  along  the 
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Earth's  magnetic  field  can  reach  an  order  of  a  hundred  kilometers  because  it  is  determined 
by  two  factors:  the  large  mean  free  path  of  the  suprathermal  electrons  and  a  small  part  of 
electron  energy  is  lost  in  one  collision.  The  extent  of  upward  elongation  will  be 
considerably  greater  than  in  the  downward  direction,  because  the  mean  free  path  above  the 
reflection  height  is  considerably  greater  than  below  it.  Consequently,  a  significant  flux  of 
suprathermal  electrons  moves  to  the  magnetosphere,  and  can  even  reach  the  magnetic 
conjugate  point. 

The  acceleration  process  due  to  the  multiple  crossing  of  acceleration  layer  is 
averaged  and  in  final  form,  depends  on  two  scalar  factors  only:  the  full  power  density  P 
absorbed  by  fast  electrons  in  the  acceleration  layer,  and  a  characteristic  parameter 
describing  the  effectiveness  of  the  acceleration  inside  the  layer  Tef.  They  are  related  to  the 
details  of  the  accelerating  process,  effective  number  of  cavitons,  their  width,  and  so  on. 

The  distribution  function  of  suprathermal  electrons  in  the  theory  of  multiple 
electron  acceleration  can  be  presented  in  the  simple  form  [Gurevich  et  al.,  1985] 

dz  "  (3) 


f0{e,z)  =  CK0  e/ 
V/ 


exp 


'  L  cos  a 


Where  Ko  is  a  modified  Bessel  function,  e  is  the  electron  energy  and  Tef  is  the 
effective  temperature  of  suprathermal  electrons,  while  a  is  the  angle  between  the  vertical 
and  geomagnetic  field.  The  normalization  constant  C  is  directly  proportional  to  the  power 
density  P  of  the  HF  wave  absorbed  by  the  suprathermal  electrons: 


C  =  - 


m 


4  ft3  T. 


(4) 


Here  8  is  the  average  fraction  of  electron  energy  lost  in  a  single  collision  with 
neutral  molecules. 

As  given  in  eq.  (3)  the  acceleration  layer  is  assumed  to  be  located  at  z  =  0,  while 
the  Lsf  factor  is  the  characteristic  relaxation  length  of  suprathermal  electrons  in  the  upward 
(+)  and  downward  directions  (-) 


Where  Nm  =  Nm±(z)  is  the  neutral  density  above  (+)  and  below  (-)  the  layer,  a,r(t) 
is  the  total  transport  cross-section  of  electron-neutral  collisions,  and  8=  ain/ <jtr,  where  a;„ 
is  the  total  cross-section  of  inelastic  collisions,  which  includes  ionization  by  electron 
impact.  In  <jtr  and  dm  the  collisions  with  all  neutral  components  are  considered. 


1.4  Comparison  of  the  Observations  with  Theory 

In  the  experiment  on  20  May  1972,  upgoing  (in  the  radar  look  direction)  fluxes  of 
suprathermal  electrons  in  the  energy  range  10  +  17  eV  at  the  heights  zj  =  256  km  and  Z2  = 
266  km  were  determined.  The  reflection  point  was  at  zo  =  285  km.  Comparison  with  the 
theory  requires  determination  of  the  electron  flux  along  radar  direction  e, : 

Jr  =  (V>er)/(V) 

We  introduce  net  upgoing  and  downgoing  fluxes  by  integrating  the  flux  Jr  over 
corresponding  angles  in  velocity  space.  Taking  the  angle  between  the  geomagnetic  field 
and  vertical  radar  ray  in  Arecibo  as  a  =  40°  we  obtain  the  upgoing  J+  and  downgoing  J 
fluxes  in  energy  interval  dft. 
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J±  =-^^0(^)fl±|(3^cosa>|exp  -  j— j4 

The  results  of  calculations  are  presented  in  Fig.  1.2  for  the  given  heights  and 
different  Tef.  It  is  apparent  from  the  figure,  that  in  the  energy  range  10  20  eV  reveal  a 

rather  flat  spectrum.  Besides,  there  is  no  strong  difference  between  the  fluxes  J  at  different 
temperatures  Tef. 

On  the  other  hand,  as  follows  from  calculations  the  fluxes  are  effectively 
diminishing  with  distance  from  the  acceleration  layer.  This  is  because  electron  energy  is 
lost  in  inelastic  collisions  with  the  neutral.  Equation  (7)  implies  that  the  fluxes  of  upgoing 
electrons  are  smaller  than  the  downgoing  fluxes  below  the  acceleration  region.  Their 
relation  depends  upon  the  angle  a  between  radar  ray  direction  and  the  geomagnetic  field  - 
for  Arecibo  at  a  =  40°  J+/J-  «  1.5.  Comparison  between  the  theory  and  observations 
presented  in  Fig.  1.3  shows  a  reasonable  agreement  between  those  two.  In  fact,  the 
behavior  of  the  spectrum  of  the  electron  flux  at  different  heights  is  consistent  with  the 
theory  for  any  of  Tef  applied.  Taking  into  account  the  absolute  values  of  upgoing  flux  J+  = 
(4-^8)xl05  el/cm2s  eV  at  characteristic  energies  IOh-15  eV  obtained  by  Carlson  et  al. 
[1982],  one  can  find  the  absorbed  power  PFs  of  the  HF  wave  converted  into  acceleration  of 
the  suprathermal  electrons 


for 
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In  these  calculations  the  Arecibo  heater  beam  was  taken  roughly  10  km  for  a  full 
absorbed  power  fVs  =  PS.  The  factor  P  was  determined  from  calculations  of  the  electron 
flux  made  at  different  distances  from  the  source,  and  from  the  results  presented  in  Fig. 1.3. 
It  is  apparent  that  the  dependence  of  Ws  on  Tef  is  not  very  strong.  Note  that  eq.  (7)  allows 
us  to  estimate  the  upgoing  electron  flux  at  z  >  500  km  which  escapes  to  the  magnetosphere 
due  to  the  absence  of  scattering  collisions. 

Figure  1.4  shows  the  height  distribution  of  plasma  line  intensity  measured  during 
the  second  experiment,  which  is  compared  with  the  height  dependence  of  the  distribution 
function  of  suprathermal  electrons  fo(e(z),z).  Here  energy  s(z)  is  determined  through  the 
measured  electron  density  distribution  (Fig.  1.1)  by  using  eq.  (2).  One  can  see  a  sufficient 
agreement  between  the  theory  and  observations.  Note,  that  the  normalized  height 
dependence  of  the  distribution  function  practically  does  not  depend  on  the  parameter  Tef 
(see  Fig.  1.4). 

1.5  Discussion  and  Conclusions 

It  was  shown  that  the  ISR  plasma  line  measurements  are  in  reasonable  agreement 
with  the  theory  of  multiple  electron  acceleration.  From  the  observational  data  and  their 
comparison  with  the  theory  it  follows  that:  In  experiments  where  powerful  radio  waves 
interact  with  the  ionosphere  a  large  number  of  suprathermal  electrons  are  generated  in  the 
energy  range  up  to  20  eV. 

The  suprathermal  electrons  are  observed  over  a  wide  altitude  region  of  the  order  of 
a  few  tens  km  both  below  and  above  the  acceleration  layer,  centered  near  the  reflection 
altitude  of  the  powerful  radio  wave.  While  not  observed,  they  must  extend  well  above  this 
altitude,  into  the  plasmasphere/magnetosphere,  and  even  to  some  degree  into  the 
magnetically  conjugate  region. 

In  power,  fVs~  5.6  *  8.5  kW  (4  -s-  6)%  of  the  full  transmitted  heater  power  goes 
into  the  acceleration  of  suprathermal  electrons.  This  energy  is  dissipated  due  to  the 
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generation  of  optical  emissions,  ionization  and  heating  of  ionized  and  neutral  components 
of  the  ionosphere  plasma.  The  characteristic  dissipation  length  for  suprathermal  electrons 
depends  strongly  on  the  altitude.  For  z  =  250  km  it  is  about  10  km,  while  for  z  ~  300  km  it 
is  about  30  km. 

The  flux  of  suprathermal  electrons  into  the  plasmasphere/magnetosphere  depends 
strongly  on  the  height  zo  of  the  acceleration  layer.  For  the  studied  cases  zo  -  285  km,  and  zo 
=  295  km,  and  Ws  =8.5  kW,  Te/=  10  eV  the  flux  into  magnetosphere  in  the  energy  range 
10  20  eV  is  about  0.9  kW.  This  estimate  is  consistent  with  earlier  calculations  of 

photoelectron  energy  and  flux  escape  into  and  loss  to  the  plasmasphere/magnetosphere 
[Mantas  et  al.,  1978]. 

We  conclude  that  plasma  line  observations  combined  with  the  theory  enables  us  to 
obtain  significant  information  about  the  acceleration  of  the  suprathermal  electrons  in 
ionosphere  modification  experiments.  For  the  first  time  the  structure  and  size  of  the 
perturbed  region  filled  with  suprathermal  electrons  along  with  the  full  power  going  into  the 
accelerated  electrons  is  determined  by  a  theoretical  calculation  consistent  with  an 
observational  data  set. 

Note  that  we  had  the  possibility  to  use  here  only  a  small  part  of  ISR  data  on 
ionosphere  existing  from  experiment  1992.  The  further  elaboration  of  this  data  and  their 
comparison  with  the  theory  is  of  a  significant  interest  and  will  be  definitely  quite  fruitful. 

We  suggest  also  these  experiments  be  repeated  to  obtain  more  information  about 
the  main  features  of  acceleration,  and  its  dependence  upon  the  reflection  height  and 
ionosphere  conditions.  Note  that  such  observations  have  as  yet  been  done  only  at  Arecibo, 
though  we  expect  that  the  latitude  dependence  of  the  effect  could  be  very  significant. 

At  the  same  time,  the  agreement  that  found  between  the  theory  and  observations 
make  it  possible  to  plan  more  detailed  future  experiments.  These  will  thus  help  to  reach  a 
much  better  understanding  of  the  physical  mechanisms  of  electron  acceleration  and 
nonlinear  processes  in  ionosphere  plasmas. 

1.6  References 

Bernhard,  P.  A.,  C.  A.  Tepley,  and  L.  M.  Duncan,  Airglow  enhancement  associated 
with  plasma  cavities  formed  during  ionosphere  heating  experiments,  J.  Geophys.  Res.,  95, 
9071,  1989. 

Carlson,  H.  C.,  V.  B.  Vickwar,  and  G.  P.  Mantas,  Observations  of  fluxes  of 
suprathermal  electrons  accelerated  by  HF  excited  instabilities,  J.  Atmosph.  Terr.  Phys.,  44, 
1089,  1982. 

Djuth,  F.T.,  J.H.  Elder,  and  K.L.  Williams,  Auroral  and  midlatitude  radar  studies: 
Radar  methods  for  improved  diagnostics  of  the  HF-  modified  and  natural  ionosphere,  Final 
Report  PL-TR-96-2070,  Geospace  Research,  Inc.,  1996. 

Dubois,  D.  F.,  A.  Hanssen,  H.  A.  Rose,  and  D.  Russell,  Space  and  time  distribution 
of  HF  excited  Langmuir  turbulence  in  the  ionosphere:  comparison  of  theory  and 
experiment,  J.  Geophys.  Res.,  98A,  17,543,  1993. 

Gurevich,  A.  V.,  Ya.  S.  Dimant,  G.  M.  Milikh,  and  V.  V.  Vaskov,  Multiple 
acceleration  of  electrons  in  the  regions  of  high-power  radio  wave  reflection  in  the 
ionosphere,  J.  Atmosph.  Terr.  Phys.,  47,  1057,  1985. 

Mantas,,  G.P.,  H.  C.  Carlson,  and  V.B.  Wichwar,  Photoelectron  flux  buildup  in  the 
plasmasphere,  J.  Geophys.  Res.,  83,  1,  1978. 

Perkins,  F.  W.,  and  E.  E.,  Salpeter,  Enhancement  of  plasma  density  fluctuations  by 
nonthermal  electrons,  Phys.  Rev.,  139,  A55,  1965. 


6 


Wang,  J.  G.,  D.  L.  Newman,  and  M.  V.  Goldman,  Vlasov  simulations  of  electron 
heating  by  Langmuir  turbulence  near  the  critical  altitude  in  the  radiation  -  modified 
ionosphere,  J.  Atmosph.  Terr.  Phys.,  59,  2461,  1997. 

Yngvesson,  K.  O.,  and  F.  W.  Perkins,  Thomson  scatter  studies  of  photoelectrons  in 
the  ionosphere  and  Landau  damping,  J.  Geophys.  Res.,  73,  97,  1968. 


Figure  Captions 

Figure  1.1  Spectra  of  the  downshifted  plasma  line  measured  at  Arecibo  13  July 
1992  at  8:30  LT. 

Figure  1.2.  Upgoing  flux  of  suprathermal  electrons  versus  their  energy  computed 
for  the  full  absorbed  power  Ws  =  10  kW  and  for  different  values  of  Tef 

Figure  1.3.  Upgoing  flux  of  suprathermal  electrons  versus  their  energy.  Solid 
curves  1,  2  and  3  correspond  to  Tef=  5,  7.5  and  10  eV  respectively,  while  the  points  with 
bars  correspond  to  the  observations  made  on  052072. 

Figure  1.4.  Plasma  line  intensities  versus  altitude.  Solid  curves  1,2  and  3 
correspond  to  Tef=  5,  7.5  and  10  eV  respectively,  while  the  points  with  bars  correspond  to 
the  observations  made  on  071392. 


Figure  1.1 
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2.  Optic  emission  from  modified  by  powerful  radio  waves  ionosphere 
Introduction 

Enhancement  of  optic  emission  was  established  in  the  first  ionosphere  modification 
experiments  (Sipler  and  Biondi,  1971;  Adeishvili  et  al.,  1976).  It  was  studied  afterwards  in 
multiple  observations  (Haslett  and  Megill,  1974;  Carlson,  1974;  Bernhardt  et  ah,  1989). 
During  the  last  years  new  interesting  results  have  been  obtained  at  HAARP  facility 
(Peterson,  1999;  Carlson,  1999),  in  sporadic  E-layer  at  Arecibo  and  at  low  duty  circle 
experiments  at  SURA  (Nazyrov  et  ah,  1999).  The  optic  emission  measurements  indicate 
directly  the  existence  of  suprathermal  electrons  accelerated  up  to  the  energies  1 0  eV  and 
more  in  modified  ionosphere. 

Other  type  of  observations,  demonstrating  the  accelerated  up  to  the  energies  20  eV 
electrons  was  performed  by  Carlson  et  ah,  1982,  Juth  et  ah,  1996,  using  incoherent 
scattering  radar  at  Arecibo.  The  suprathermal  electrons  were  seen  in  these  experiments  in  a 
wide  region  of  the  order  100-km  around  the  reflection  point  of  the  heater  wave.  Some 
indications  of  the  significant  effect  of  suprathermal  electrons  on  SEE  emission  exist  as 
well  (Frolov,  1998). 

Theoretical  considerations  from  the  very  beginning  connected  the  enhancement  of 
optic  emission  with  the  acceleration  of  electrons  in  the  Langmuir  plasma  turbulence  layer 
near  the  reflection  point  of  powerful  HF  O-wave.  Though,  it  was  indicated  as  well,  that  the 
mostly  intensive  red-line  emission  (630  nm)  partly  could  be  determined  by  the  direct 
effective  heating  of  the  main  bulk  of  electrons  (Carlson,  Mantas,  1996),  or  it's  significant 
part,  captured  inside  nonlinear  density  depletions  -  "striations"  generated  by  the  upper- 
hybrid  turbulence  (Gurevich  and  Milikh,  1998).  The  theory  of  Langmuir  plasma 
turbulence  generated  near  the  reflection  region  of  pump  O-wave  connects  the  electron 
acceleration  process  with  the  modulation  instability  and  with  the  creation  of  nonlinear 
Langmuir  cavitons  -  density  depletions  filled  with  Lagmuir  plasma  oscillations  (Wong  and 
Stenzel,  1975).  The  fast  electrons,  passing  trough  the  turbulent  region  filled  with  cavitons 
obtain  the  energy  from  the  trapped  caviton  plasma  oscillations  in  a  diffusive  way 
(Eastebrook  et  al.,  1975;  Morates  and  Lee,  1977;  Andreev  et  al.,  1980;  Wang,  Goldman 
and  Newman,  1 997).  The  main  feature  of  this  mechanism  is  that  the  energy  is  gained  only 
by  fast  electrons  whose  velocity  is  high  enough  v  >  a/co,  where  a  is  a  scale  of  a  caviton,  co 
is  the  heater  wave  plasma  frequency.  Low  energy  electrons  oscillate  in  the  caviton 
adiabatically  and  do  not  get  any  additional  energy.  It  should  be  emphasized,  that  the 
collapse  of  cavitons  does  not  effect  strongly  the  acceleration  process  of  fast  electrons 
(Wang  et  al.,  1997).  In  the  same  time  collapse  can  input  a  significant  part  of  wave  energy 
to  the  thermal  component  (Du  Bois  et  al.  1993,  Hanssen  et  al.  1997).  To  calculate  the  flux 
of  energy,  which  goes  to  fast  electrons  the  detailed,  kinetic  theory  is  needed,  which  takes 
into  account  electron  collisions. 

For  example,  the  number  and  energy  distribution  of  fast  electrons  in  ionosphere 
conditions  is  affected  by  the  so-called  "multiple  acceleration"  (Gurevich  et  al.,  1985, 
Vaskov  et  al.,  1983).  The  multiple  acceleration  process  is  determined  by  the  fact,  that  the 
Langmuir  acceleration  take  place  in  a  strongly  excited  layer,  embedded  in  a  weakly 
disturbed  ionosphere  plasma.  After  the  fast  electrons  leave  the  acceleration  layer  it  collides 
with  a  neutral  molecules  in  nondistrubed  plasma.  Due  to  collisions  some  electrons  could 
return  to  the  acceleration  layer  and  get  there  the  energy  once  more.  The  process  could  be 
repeated  many  times,  thus  creating  a  "multiple  acceleration".  Directly  this  process 
determines  the  high-energy  tail  of  electron  distribution  function  and  a  wide  spread  of  fast 
electrons  around  the  acceleration  region. 
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Though,  understanding  of  the  main  mechanisms  of  the  acceleration  process  in 
Langmuir  turbulence  seems  significantly  clear  enough,  both  the  theory  and  observations  in 
ionosphere  modification  experiments  up  to  now  were  not  developed  to  the  state  when  a 
detailed  quantitative  comparison  between  them  could  take  place.  That  means,  that  some 
fundamental  problems  of  modification  experiments  remain  non-answered.  One  of  them  is 
the  intensity  and  space  distribution  of  the  observed  optic  emission. 

In  the  present  work  we  intend  to  develop  the  existing  theory  and  apply  it  to  obtain 
the  detailed  comparison  with  the  main  part  of  the  existing  experimental  data.  The  paper  is 
constructed  as  follows.  In  section  2  heating  of  electrons  by  both  in  upper  hybrid  and  in 
Langmiur  turbulence  is  considered.  It  is  demonstrated  that  the  calculations  of  electron 
acceleration  in  Langmiur  resonance  layer  without  electron  collisions  lead  to  controversial 
result.  In  section  3  the  theory  of  multiple  acceleration  is  generalized.  The  energy  losses  of 
electrons  on  the  excitation  of  molecular  vibration  level  are  taken  into  account.  In  section  4 
the  model  calculations  describing  the  general  peculiarities  of  the  suprathermal  electron 
space  energy  distribution  are  performed.  The  section  5  describes  the  detailed  comparison 
of  the  theory  with  the  main  observational  data  of  the  artificial  optic  emission  in  F-layer. 
The  strong  enhancement  of  the  optic  emission  in  sporadic  E-layer  modification  is 
discussed  in  sections  6.  In  conclusion  we  formulate  the  nowadays  state  of  the  problem. 
The  strongly  inhomogeneous  character  of  the  disturbed  region  is  stressed. 


1 .2  Heating  of  electrons  in  the  resonance  layer 

In  the  vicinity  of  reflection  point  of  powerful  O-wave  in  ionosphere  both  heating 
and  acceleration  of  plasma  electrons  takes  place. 

As  is  well  known  for  vertically  propagating  pump  O-wave  the  resonance  region  is 
situated  at  the  heights  Z  between  upper  hybrid  ZUb  and  Langmuir  Zl  resonances 
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In  this  region  the  frequency  of  the  pump  wave  co  coincide  with  the  frequency  of  the  natural 
oscillations  in  ionosphere  plasma 
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Here  N(Zl )  -  electron  density  at  the  Langmuir  resonance,  where  the  Langmiur  turbulence 
take  place.  And  ZUh  is  defined  by  resonance  condition  for  upper  hybrid  resonance 
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In  the  vicinity  of  Zuh  the  upper  hybrid  (UH)  turbulence  is  effectively  developing.  The 
width  of  resonance  region  AZ  depends  on  the  pump  wave  frequency  and  plasma  density 
gradient  L : 


A  Z  =  ZL-Zuh 
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coc  is  electron  gyroferequency.  For  the  ionosphere  F-layer  conditions  usually  AZ  ~  2  +  5 
km. 

Conventional  collision  absorption  of  radio  waves  in  F-region  is  very  weak  -  less 
than  1  dB  [Ginzburg,  1967].  If  plasma  turbulence  in  the  resonance  region  is  excited  the 
absorption  become  abnormally  strong  -  up  to  10+20  dB  [Gurevich,  1978]. 

Note,  that  the  Langmuir  turbulent  layer  (LT)  is  placed  in  the  upper  part  of 
resonance  region  close  to  the  reflection  point  (2),  while  close  to  the  upper  hybrid 
resonance  (3)  lays  upper  hybrid  turbulent  layer  (LH).  Acceleration  of  electrons  take  place 
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in  Langmuir  resonance  layer  (we  do  not  speak  here  about  singled  out  case  of  multiple 
gyroresonance  frequency  co  =  ncoc).  The  UH  turbulence  lead  to  formation  of  striations.  An 
effective  anomalous  absorption  of  the  pump  wave  take  place  in  this  region,  which  lies 
below  the  Langmuir  turbulence  layer.  It  means  that  the  UH  anomalous  absorption 
(absorption  on  striations)  can  affect  significantly  the  heating  and  acceleration  processes  in 
the  Langmuir  turbulence  layer  while  diminishing  the  pump-wave  power. 


Heating  of  electrons  in  UH  turbulence  region 

Anomalous  absorption  is  determined  by  the  excitation  of  UH  waves  due  to  direct 
transformation  of  the  pump  wave  on  the  density  gradient  in  striations.  UH  waves  trapped 
in  striations  heat  electrons  and  thus  determine  the  depth  form  and  number  density  of 
striations.  Nonlinear  theory  of  the  electron  heating  and  anomalous  absorption  on  striations 
was  developed  by  Gurevich  et  al.  1996.  The  results  of  the  theory  are  in  a  rather  good 
agreement  with  observations. 

Following  to  the  theory  of  UH  anomalous  absorption  we  determined  the  simple 
relation  for  the  losses  of  pump  wave  power  P  in  the  UH  turbulence  region  on  its  way  to 
Langmuir  turbulent  layer. 
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Here  K  is  an  attenuation  coefficient  in  dB.  Coefficient  Q  is  proportional  to  the 
relative  depth  of  striations.  It  depends  on  pump  wave  power,  frequency,  plasma 
parameters.  The  state  of  the  striation  development  is  also  significant.  For  the  well 
developed  striations  Q  ~  1 . 

1 .3  Heating  of  electrons  in  Langmuir  turbulent  layer 

The  excitation  of  Langmuir  turbulence,  it  steady  state  and  caviton  formation  in  the 
field  of  powerful  radio  wave  in  conditions  close  to  ionosphere  was  studied  numerically. 
We  will  discuss  the  result  of  the  calculations  of  Hanssen  et  al  and  Wang  et  al  in  order  to 
determine  the  thermal  heating  of  electron  by  Langmuir  turbulence.  The  parameters  used  in 
calculations  are  presented  in  Table  1 . 


Hanssen  et  al 
(1992) 

Wang  et  al  (1997) 

Pump  wave  frequency  (MHz) 

f=  4.01 

mmmm 

Pump  wave  field  Erms  (V/m) 

1.26 

1.6 

Plasma  density  no  (cm'J) 

Electron  temperature  (eV) 

E] 

0.1 

Pump  wave  energy  Wrms=Erms2/4m0To 

fflHH 

Hanssen  et  al.  (1992)  considered  one-dimensional  Zakharov  system  of  equations 
describing  nonlinear  development  of  plasma  density  and  electric  field  fluctuations  under 
the  action  of  homogeneously  oscillating  electric  field  of  the  pump  wave.  Electron  and  ion 
collisions  and  Landau  damping  were  taken  phenomenologycally,  using  terms  from  linear 
theory.  The  established  average  energy  of  longitudinal  plasma  waves,  concentrated  mostly 
in  cavitons  grew  up  70  times  in  comparison  with  the  pump  wave  energy.  It  gives 
W^rms  x  0.19. 

The  power  density  dissipating  by  plasma  oscillations  to  the  heating  of  electrons  is 
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Where  s  is  coordinate  along  the  Earth  magnetic  field,  ve  -  electron  collision 
frequency,  co  -  plasma  wave  frequency  and  Erms  =  4moTe  W^rms-  Following  Hanssen  et  al. 
let  us  suppose  that  the  whole  width  of  the  turbulent  layer  is  of  the  order  of  pump  wave 
length  near  reflecting  point  L  ~  A,  «  200  m.  We  obtain  that  the  full  pump  wave  power 
dissipating  in  Langmuir  cavitons  to  the  heating  of  plasma  electrons: 
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That  value  is  quite  significant.  Really,  if  we  neglect  anomalous  absorption  in 
striations  below  turbulent  layer,  then  to  create  the  electric  field  Erms  =  1 .26  V/m  (see  Table 
1)  at  the  main  Eiry  maximum  near  reflection  point,  the  pump  wave  power 
60(sina)4/3(100km/L)l/3/jKW/m2  is  needed.  Here  a  is  the  inclination  angle  of  Earth 
magnetic  field.  We  see  that  at  a  ~  1  about  5  dB  of  pump  wave  is  absorbed  due  to.  the 
thermal  heating  of  electrons  in  the  cavitons.  Note  that  the  calculated  pump  wave 


absorption  is  effectively  growing  with  latitude  (a'4/3).  So,  more  than  60%  of  absorbed 
pump-wave  energy  in  caviton  turbulence  is  going  directly  to  the  thermal  component  of 
electrons. 

Wang  et  al.  1997  elaborated  one-dimensional  kinetic  model.  Not  only  field 
equation  but  also  Vlasov  kinetic  equation  for  noncolliding  electrons  was  integrated 
simultaneously  in  the  model.  The  Langmuir  waves  instability  was  excited  and  established. 
Energy  level  of  average  plasma  oscillations  grew  up  about  30  times  in  comparison  with 
pump  wave  field  reaching  Wnns  =  0.14.  We  see  that  in  spite  of  quite  different  approach,  the 
numerical  value  of  saturated  energy  of  plasma  oscillations  is  close  enough  in  both 
calculations. 

The  power  density  dissipating  to  the  thermal  component  of  electrons  is  again 
determined  by  expression  (6).  The  width  of  turbulent  layer  was  supposed  to  be  Xm  ~  180 
m.  From  (6)  we  again  obtain  the  heating  power  Wj=  40  p.  kW/m2,  what  means  that  about 
2.5  dB  of  pump  wave  energy  goes  to  the  heating  in  cavitons. 

We  see  that  when  Langmuir  turbulence  is  effectively  excited  the  pump  wave 
energy  losses  to  the  heating  of  electrons  is  quite  significant,  compatible  (though  remaining 
less)  to  the  anomalous  absorption  in  UH  layer. 

Significant  advantage  of  Wang  et  al.  1997  kinetic  model  is  the  possibility  to  see 
directly  not  only  saturation  of  plasma  oscillations,  but  acceleration  of  suprathermal 
electrons  as  well.  The  result  of  calculations  shows  that  a  tail  of  initial  Maxwelian 
distribution  function  is  growing  in  time.  An  established  suprathermal  stationary  tail  is 


/,  =8.3x10" 


32 V* 

\  v  ) 


22L  4<  — <25 

Vt.  VTe 


(7) 


shown  to  be  approximated  by  power  -  law: 


Here  VTe=(Te/m),/2  _  thermal  velocity.  The  one  dimensional  maxwellian  part  of 


distribution  function  has  standard  view: 
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The  full  number  of  suprathermal  electrons  is  not  large  about  2%  of  the  thermal  electrons. 


13 


Let  us  determine  now  the  transport  of  energy  by  fast  electrons.  Integrating 
distribution  function  (7)  we  can  calculate  the  full  power  transported  out  of  the  turbulent 
layer  by  accelerated  suprathermal  electrons 

h\  ukw 

Wa=  2  [  —  mv2vfsdv  «  2.68n0TcVrc  « 1700^-— —  (9) 

41 't,  2  m 

We  see  that  the  power  transported  by  suprathermal  electrons  is  20  times  larger, 
than  the  full  power  of  the  pump  wave  ( Wp  »  96jukW/m2).  This  result  is  obviously 
controversial.  It  means,  that  in  spite  of  the  fact  that  mean  free  path  of  the  suprathermal 
electrons  in  ionosphere  F-layer  is  larger  than  the  width  of  the  turbulent  region  Am  «  1 80  m, 
the  electron  collisions  could  not  be  neglected.  The  flux  of  the  energy  out  of  the  layer  is 
determined  by  collisions  _  suprathermal  electrons  not  only  leave  the  layer,  but  due  to 
collisions  they  can  return  back.  Thus  only  the  full  balance  between  accelerated  inside  the 
layer  and  returning  backward  suprathermal  electrons  can  determine  correct  relation 
between  the  pump  wave  energy  losses  and  acceleration.  Just  this  balance  is  studied  in  the 
theory  of  multiple  acceleration.  This  theory  in  modernized  form  would  be  considered 
below. 

Note,  that  the  controversial  result  of  electron  acceleration  in  non  collisional 
kinetics  shows,  that  the  theory  of  plasma  turbulence  saturation  should  be  revised  as  well: 
electron  collisions  outside  the  Langmuir  turbulent  layer  should  be  taken  into  account  to 
organize  a  correct  flow  of  plasma  energy.  In  other  words  boundary  conditions,  which 
connect  turbulent  and  nonturbulent  regions,  could  affect  the  process  of  Langmuir  waves 
saturation  inside  the  turbulent  layer. 


1 .4  Kinetic  theory  of  multiple  acceleration 

To  find  the  distribution  function  we  have  to  solve  kinetic  equation  inside  and 
outside  acceleration  layer  and  match  solutions. 

Inside  the  layer  kinetic  equation  has  a  form: 


df  df  e  df 
—  +  v„— - E-^—  =  0 


dt 


dz 


m 


dv„ 


Here  E  is  electric  field  of  Langmuir  cavitons  excited  by  powerful  radio  wave 

E(z,t)  =  ^E{z)(e-',°'  +c.c) 

The  solution  of  the  equation  (10)  has  the  form: 
f(z,  spt)  =  F(z  - Z(t),el ,  -  s(t)) 

Here  eft),  Z(t)  the  trajectory  of  a  fast  particle  in  electric  field  E: 


(10) 

GO 

(12) 


—  =  V,  —  =  —E,  V  =  j2e/m 
dt  dt  m 

To  solve  our  problem  we  need  to  determine  the  changing  of  distribution  function 
Af  when  a  fast  particle  crossing  the  acceleration  layer.  Averaging  relation  (12)  on  fast 
oscillations  we  find  the  changing  of  distribution  function  after  passing  through  the 
acceleration  layer). 


A/: 


\__d_ 

2  d£n 


(13) 


According  to  (1 1)  after  averaging  on  fast  oscillations  in  a  fix  spatial  point  energy  gain  is 
given  by  a  formulae: 
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A£'(z,r)=-A£,COS 


cot - +  y/ 

V  vll  J 


(14) 


+co 

A£{e^)-e\Ek\,  Ek  =  §  E{z)e~'k:  dz, 


i  03 

k-  —  =  co 


i  m 

2a, 


For  example,  for  Gauss  like  caviton  E(z)=  Eoexp(-(z/a)2)  we  have  from  (11): 


As  =  7c^eE0a  exp(-  ka 2 ) 


(15) 


We  see  from  (15)  that  only  electrons  with  high  enough  energy 


£j,  >  s.=  m(coa)2  =  Te 


f a > 

2 

c  \ 

a 

U  ; 

IaJ 

»1 


(16) 


where  De  is  a  Debay  length,  could  be  accelerated.  It  is  easy  to  see  that  relation  (16)  has  a 
universal  character  and  does  not  depends  on  a  real  form  of  a  caviton.  It  is  obvious  that  if 
we  have  some  number  of  noncorrelated  cavitons  the  average  energy  changing  is  equal  to  a 
sum 

(Af2)=Z(A*2) 

i 

Now  let  us  consider  the  solution  outside  acceleration  layer  at  z>0  and  z<0.  Kinetic 
equation  has  a  form: 

(17) 


df  df  _ 
—  +  v«—  =  -S 
dt  dz 


Here  f(t,v,[c,z)  is  the  distribution  function  of  electrons,  S  -  collision  integral,  p  _ 
cosine  of  an  angle  between  the  velocity  of  electron  and  z  direction.  The  collision  integral 
consists  of  two  parts:  inelastic  So  and  elastic  Si.  The  collision  integral  depends  strongly  on 
energy.  For  energies  e<s*  (e*«4  eV)  the  main  collision  crossection  is  defined  by 
oscillational  levels  and  inelastic  collision  integral  has  a  diffusive  form.  But  for  energies 
s>e*  the  collision  crossection  is  determined  by  optical  levels.  Thus  the  collision  integral 
could  be  presented  in  a  form: 


f 


S0(F0)  =  -1r 

os 


,df 


R^- 
V  d£ 


S0(F0)  =  v0(s)F0 


£  <£ 


v0  =  va0Nm 


£>  £ 


(18) 


v0  =vv,Nm 


Si(f)  =  Vl(£)f], 

a,  =ate{£-£)+YuQ(£-£') 

t 

The  <7o  isthe  sum  of  effective  crossection  of  inelastic  collisions  with  neutral  atoms  N„,  -  the 
total  density  of  neutrals.  at  is  transport  electron  crossection  of  electron  elastic  collisions 


cr 


0  — 


N. 


mk 


N„ 


(19) 


The  averaged  part  of  total  energy  of  an  electron  loosing  in  one  collision  is  small;  it  is 
defined  by  parameter 


8  =  vjv ,  =<7,1(7, 


(20) 
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High  electron  energies  (e  >  £*) 

Because  of  (20)  we  can  expand  the  equation  (17)  on  small  parameter  5.  At  first 
approximation  the  distribution  function  does  not  depend  on  p.  Taking  into  account  first 
approximation  one  can  find: 


dt  3  dz 


(21) 


dFr 


dz 

Boring  in  mined  the  boundary  conditions  (13)  at  z  =  0  in  stationary  conditions  from  (21) 
one  can  find: 


dF0  3v,  d 

— -  =  -sign(z ) — - - 

dz  9vs  d£ 


j(Af  (^||  ))d£l} 


(22) 


It  is  easy  to  see  that  for  collision  integrals  in  the  form  (18)  the  equation  (22)  depends  on 
energy  parametrically  only.  So  we  can  rewrite  equation  (18)  in  the  region  of  energies 


d2F0 

d£2 


=  3  SF0, 


£  =  [^dz=l°xNm(z)dz 


(23) 


Taking  into  account  that  5  is  a  slow  changing  function  it  is  possible  to  solve  (23)  in  WKB 
approximation 


F0(£,z)  =  F0(£)\ 


"A(0)" 

K  1 

>  dz 

Az). 

expj  - 

!laz) 

(24) 


/,  (A,  2)  =  sign(z)(3S)i//2  F0  (e,  z ) 

Here  Le  is  a  relaxation  scale  of  fast  electrons  with  given  energy 


Middle  electron  energies  (Te<£<£*) 

In  the  region  of  energies  £<  £*  the  kinetic  equation  (21)  takes  the  form: 


d2Fn 


d ^  d£ 


d  ( Rdpo 


(25) 


v  de  j 

Since  the  solution  of  kinetic  equation  in  the  region  £  >  e*  has  a  form  (24),  to  match  this 
solution  with  solution  in  region  e<  £*  we  had  to  search  the  solution  for  £  <  £*  in  a  form 


=g(A)exp(- 


L 

zr  dz 

1 

]i,w  J 

The  equation  (25)  takes  the  form: 


d  ^ rQL' 

9e 


v 


8e 


Taking  into  account  that 


v0v,  «  Rv{ 


One  can  see  from  (26)  that  at  first  approximation 
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£  <  £ 


F0  (£•)  =  F0  (£■* )  =  const  £<£*  (27) 

The  distribution  function  Fo(e)=Fo(s,z)|z=o  according  to  boundary  condition  (13)  is 
determined  by  diffusion  equation: 


1  d  f£^.%=/0-/(‘V), 


lz  dsY  'eJJd£ 


8/2e 

and  Fq(e)  =  Fa{£)  =  const,  £<£ 

We  introduce  here  effective  temperature : 


Tcff  {£)  =  kAe,  As  =  |  -  |(Af 2  («■,  ))det  , 


F0=f0-f(0)(£ ) 


3A  -y 

K  —  ——^=8  /4 


2V2 


(28) 


(29) 


here  Af  is  mean  electron  energy  increasing  in  accelerating  layer.The  effective  temperature 
increase  gradually  with  increasing  energy  and  saturate  at  the  level. 

Tc/rm=K(n/2)//2Asm  £m<£l 

Flere  n  is  the  number  of  cavitons  in  accelerating  layer. 

The  intensity  of  electron  acceleration  is  defined  by  parameter 

'dfm{£) 


T 

.  _  Fff 


T  =  -fm(£ ) 


T  '  7  ds 

For  weak  acceleration  (y<l)  the  distribution  function  of  electrons  does  not  change 
substantially 


F0(£)  = 


f  T  \2 

V  T  , 


f(0)(£) 


In  opposite  case  of  strong  acceleration  y  »  1  the  equation  (29)  in  WKB  approximation 
has  the  following  solution: 

f0(£)  =  fm(£)  £<£ih 

fo(£)  =  M£,h)<  ^  expi-  | - |  £>£„ 


Yi ,  ,  -\/ 

__ ^ 

s- 

2£Teffm(£)8'A(s) 

exp 

[  ds  j  e  P 

1 

! 

C 

'ST 

1 _ 

(30) 


The  value  sth  in  (30)  is  defined  by  matching  of  distribution  functions  in  saddle  point 

Fej]m  —F(£lh). 

In  the  case  of  strong  acceleration  in  the  limit  Teff  *  const,  8  «  const  the  distribution 

function  (30)  has  the  form: 


m=cK0 


c  \ 
£ 


\T‘ff  J 


(31) 


Flere  Ko( x)  is  a  modified  Bessel  function.  The  number  of  fast  electrons  is  defined  by  the 
integral 

n!  =  jf(£)d\ 

But  the  normalized  constant  C  is  determined  by  the  microwave  power  absorbed  by  fast 
electrons.  The  absorbed  power  consists  of  two  terms 
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The  first  term  Pi  in  (32)  is  responsible  for  loses  on  rotational  levels  of  molecules  and  the 
second  one  P2  describes  the  electron  acceleration. 


The  work  will  be  prolonged  as  it  is  formulated  in  the  Introduction.  Its  final  text  will 
be  presented  in  the  Final  Report. 
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3.  Drift  effects  on  the  scale  of  field-aligned  irregularities 


Irregularities  in  the  F-region  of  the  ionosphere  are  strongly  elongated  along  the 
magnetic  field  due  to  dramatic  difference  in  transport  coefficients  along  the  magnetic  field 
and  transverse  to  it.  Their  scales  along  the  magnetic  field  L||  in  stationary  state  were 
determined  previously  neglecting  plasma  drift  effects  both  in  linear  [Gurevich,  1978  and 
many  other  publications]  and  nonlinear  [Gurevich  et  al.,  1995]  approximations.  One  can 
suppose  that  drifts  could  substantially  diminish  L||  for  irregularities  of  sufficiently  small 
transverse  (to  the  magnetic  field)  scales  in  the  drift  direction. 

Transport  equations  along  the  magnetic  field  can  be  written  in  the  form  [Gurevich, 


1978] 

dT 

S  | 

( 

dT 

-V — 

H - 

Kw 

dx 

dz' 

ii 

V 

dz 

-SvT  +  Q  =  0 


(1) 


dx  dz 


( 


a 


dN\  d 


dz 


dz 


D, 


rdJ 
dz , 


-yN  =  0 


(2) 


Here  T  and  N  are  disturbances  of  temperature  and  density  in  irregularities,  D,  D  ,  and  k 
with  index  ||  are  diffusion,  thermal  diffusion,  and  heat  conductivity  along  the  magnetic 
field  respectively,  y  is  the  recombination  coefficient,  8  is  the  average  part  of  electron 
energy  loss  under  collisions  with  ions  and  neutrals,  and  \nu  is  the  electron  collision 
frequency.  In  (1),(2)  we  supposed  that  irregularities  are  induced  by  heating  source  Q,  we 
neglected  all  transport  across  the  magnetic  field  (supposing  transverse  scale  to  be  not  too 
small)  and  took  into  account  drift  across  the  magnetic  field  with  the  velocity  V. 

Without  drifts  longitudinal  scale  of  irregularities  is  determined  by  the  longitudinal 
diffusion  length  LN=(D\/y)l/2  and  the  longitudinal  heat  conductivity  length  Z.7=(Kj|/5v)1/2 
[Gurevich,  1978].  In  the  F-region  of  the  ionosphere  Lt  »  Ln  [Gurevich,  1978],  and 
longitudinal  scale  of  irregularities  is  determined  mainly  by  the  heat  conductivity  length  Lt- 
In  case  of  the  heating  source  Q  produced  by  the  powerful  HF  radio  waves  its  scale  along 
the  magnetic  field  is  much  smaller  LT ,  and  thus  its  heating  can  be  treated  using  Dirac  delta 
function  8(z). 

Let  us  consider  heat  transfer  equation  (1).  Introducing  undimensional  variables 

T  ,  x  z  V  Q 

?  =  —  >  v  =  —  >  9  = 


T  =  ■ 


^  T  ’ 
*~*V 


Lr 


K 


Tr.^KV 


where  Lv  =  Vo/Sv,  T0  is  unperturbed  electron  temperature,  and  Vo  is  a  characteristic  drift 
velocity,  we  can  rewrite  equation  (1)  in  the  following  form 
dz  S2 t  .  n 

~v^7  +  TT-r  +  ^te)  =  0  (3) 

dd,  dq 

In  case  of  constant  drift  velocity  v=l  one  can  easily  found  solution  of  (3)  for  harmonic  in  £, 
irregularities  x  cc  exp(ik^)  as 

r* =  2==rexp(-  ^x+iks)  (4) 

2vl  +  ik 

where  qk  is  the  respective  harmonic  component  of  the  heat  source.  From  this  solution  one 
can  find  that  irregularities  with  transversal  scale  L  <  Lv  in  the  direction  of  plasma  drift 
have  longitudinal  scales  decreasing  as  (L/Lv)1/2  as  well  as  the  efficiency  of  their  excitation. 
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In  the  F-region  characteristics  scale  Ly  is  of  order  of  kilometer  for  the  drift  velocity  of 
order  of  few  ten  meters  per  second. 

Above  conclusion  is  valid  only  if  heating  source  has  a  fixed  position  in  the  ionosphere 
(say,  forming  by  the  radio  wave  beam).  From  the  other  hand  small  scale  (few  meters) 
irregularities  -  striations  [Gurevich  et  al.,  1995]  are  excited  inside  the  beam,  and  their 
source  can  be  considered  as  drifting  with  plasma.  In  this  case  gradients  of  the  drift  velocity 
becomes  essential.  Supposing  linear  gradient  of  drift  velocity  along  the  magnetic  field  we 
can  rewrite  equation  (1)  using  undimentional  variables  as 


3r  Qz  d2z 

H - r~z  +  qS(g)  =  0 


where  characteristic  scale  Ly  is  determined  now  as 


(5) 


L  _  Lr  dVo 
v  Sv  dz 


Vanishing  at  infinity  solution  of  (5)  for  harmonic  component  t  oc  exp(ik^)  is 


(ik)  Aqk 


k~  f  V  i\ 
,  k) 


2  Ai'\ 


(6) 


Analizing  behavior  of  solution  (6)  for  different  values  of  k  one  can  find  that  longitudinal 
scale  remains  Lt  for  transversal  scales  of  irregularities  L  >  Ik  decreasing  for  smaller  scales 
approximately  as  (2lg(Lv/L))'' .  Efficiency  of  excitation  also  decreases  for  small  scales  as  a 
factor  of  k~u . 

In  the  F-region  characteristics  scale  Ly  is  of  order  of  10  meters  for  drift  velocity  gradient 
dVo/dz  of  the  order  of  several  cm/s  per  kilometer.  Thus,  we  see  that  drift  gradients  could 
essentially  affect  stationary  state  of  striations  of  meter  transverse  scales 

In  our  future  work  we  intend  to  take  into  account  effects  of  the  drift  on  the 
nonlinear  structure  of  striations  both  already  existing  in  ionosphere  and  induced  by  the 
action  of  powerful  radio  wave.  The  final  results  would  be  presented  in  our  Final  Report. 
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